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Technical  Report  Summary 


This  report  is  divided  into  four  sections.  Section  1  concerns  inferences  on 
the  crust  and  upper  mantle  structure  of  the  Central  Nevada  Test  Site  from 
teleseismic  P-wave  residuals  and  its  implications  for  the  Faultless  magnitude 
versus  yield  anomaly.  Hot  Creek  Valley,  the  site  of  CNTS,  is  dominated  by  a 
series  of  nested  calderas  which  range  in  age  from  25  to  35  m.y.,  10  to  20  m.y. 
older  than  the  Silent  Canyon  caldera  at  Pahute  Mesa.  P-wave  residuals  relative 
to  a  Great  Basin  average  show  that  Hot  Creek  Valley  is  relatively  fast,  suggesting 
that,  as  for  the  case  of  Silent  Canyon  at  NTS,  a  high  velocity  vestage  of  a  root  to 
the  Hot  Creek  calderas  exist  at  depth.  Comparing  travel-time  residuals  of  explo¬ 
sions  in  Hot  Creek  Valley  and  Silent  Canyon  relative  to  an  explosion  at  Yucca 
Flat  indicates  slower  relative  velocities  beneath  CNTS  than  at  NTS.  An  explana¬ 
tion  for  this  is  that  the  body  under  CNTS  is  probably  older  and  thus  has  had 
more  time  to  be  assimilated  into  the  surrounding  mantle.  Ray  tracing  through 
the  proposed  mantle  structure  beneath  Pahute  Mesa  (Spence,  1974)  results  in  a 
shadow  zone  for  European  stations.  A  large  percentage  of  the  stations  reporting 
magnitudes  for  the  Nevada  explosions  are  from  Europe.  The  effect  is  much 
smaller  at  Hot  Creek  Valley  suggesting  that  the  magnitude  vs  yield  bais  might 
result  from  using  magnitude-yield  curves  for  the  larger  explosions  based  on  data 
from  Pahute  Mesa.  We  have  now  obtained  the  WWSSN  data  for  FAULTLESS,  GREE¬ 
LEY,  and  DUMONT  and  are  plotting  amplitude  vs  distance  curves  to  see  if  in  fact 
the  expected  caustic  and  shadow  zone  predicted  exist  in  the  data. 

Section  2  is  a  manuscript  which  is  in  press  in  the  Journal  of  Geophysical 
Research.  In  this  study,  spectra  for  24  Mammoth  Lakes  earthquakes  greater 
than  1  x  1021  dyne-cm,  including  the  Mi  6  event  occurring  at  1450  UT  on  May  27, 
1980,  have  been  determined.  The  short-period  portion  of  the  spectra  (0.5  -  10.0 
Hz)  have  been  determined  from  spectral  analysis  of  the  locally  recorded  strong- 
motion  and  broad-band  digital  data.  A  relationship  has  been  derived  between 
the  Marshall-Basham  surface-wave  magnitude  and  seismic  moment  allowing  us 
to  determine  the  seismic  moment  directly  from  time  domain  amplitude  meas¬ 
urements  of  the  surface  waves  recorded  at  regional  stations.  This  gives  an  esti¬ 
mate  of  the  seismic  moment  independent  of  that  obtained  from  the  S-wave  spec¬ 
tra  and  one  based  on  lower  frequencies  (0.05  -  0.125  Hz).  A  comparison  of  the 
seismic  moments  indicate  that  for  any  one  event,  no  significant  difference  in  the 
spectral  level  exists  between  the  values  determined  at  surface  wave  frequencies 
(0.05  to  0.125  Hz  )  and  the  S-wave  corner  frequency  (  1  to  5  Hz  ).  The  events  stu¬ 
died,  and  in  particular  the  5/27/80--1450  Mi  8  event,  do  not  show  peaked  spec¬ 
tra  predicted  from  theoretical  analysis  of  the  ground  motion  resulting  from  a 
fluid-driven,  propagating  tensile  crack.  We  have  found  nothing  in  the  spectra  of 
these  events  to  distinguish  them  from  "tectonic"  earthquakes.  The  absence  of  a 
peaked  spectrum  does  not  rule  out  dike  injection  as  the  source  of  these  earth¬ 
quakes,  but  it  does  rule  out  a  particular  model  of  magmatic  intrusion,  namely 
that  of  Chouet  (1981)  and  Aki  (1984),  for  a  low  viscosity  fluid. 

Section  3  is  a  draft  of  a  manuscript  prepared  for  submission  to  the  Geophy¬ 
sical  Journal  of  the  Royal  Astronomical  Society.  In  an  earlier  study,  Wyss  and 
Brune  (1971)  compared  short  period  and  long  period  excitation  of  events  in  the 
California  region.  They  found  that  events  occurring  to  the  north  of  Bishop, 
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California  indicated  higher  apparent  stress  than  events  along  the  San  Andreas 
fault.  In  this  study  we  have  looked  in  detail  at  the  relative  excitation  of  short 
and  long  period  energy  in  this  region,  from  aftershocks  of  the  Mammoth  Lakes 
sequence.  We  have  compared  the  body  wave  magnitude,  m6 ,  and  surface  wave 
magnitude,  Ms,  and  find  that  Mammoth  Lakes  events  discriminate  from  nuclear 
explosions  occurring  at  the  Nevada  Test  Site  approximately  150  km  to  the 
southeast.  We  have  used  local  strong-motion  and  broad-band  digital  recordings 
to  compute  the  source  spectra  of  37  earthquakes.  For  events  in  common  with 
Archuleta  et  al.  (1982)  our  moments  are  in  agreement.  We  have  developed  a 
moment-surface  wave  magnitude  relationship  which  allows  us  to  rapidly  deter¬ 
mine  the  seismic  moment  at  surface  wave  periods  from  time  domain  amplitude 
measurements  on  the  long-period  seismogram.  These  determinations  of 
moment  in  almost  all  cases  agree  closely  with  the  moments  determined  from 
the  spectra.  Similar  studies  in  southern  California  show  long-period  moments 
exceeding  short-period  moments  by  a  factor  of  4  to  10.  For  six  events  wp  have 
compared  the  magnitude  estimates  mj,  Mi,  and  Ms  with  the  spectral  amplitude 
at  their  appropriate  frequencies.  Of  the  four  shallow  events  compared,  three, 
those  at  801471857,  801481901,  and  801490516,  are  consistent  with  the  cj~2 
model.  The  stress  drops  of  these  events  are  nearly  similar  but  the  source 
dimension  of  the  801471857  event  is  nearly  twice  that  of  the  other  two  events, 
which  explains  its  larger  moment  and  greater  surface  wave  excitation.  The 
fourth  shallow  event  had  an  anomously  low  surface  wave  excitation  compared  to 
other  events  and  a  lower  stress  drop.  The  remaining  two  events,  at  depths  of  8 
and  12  kilometers,  have  similar  stress  drops  sued  source  dimensions  to  the 
801471857  event.  However,  the  8  km  deep  event  has  anomously  large  surface 
wave  excitation,  and  the  mb  of  the  deepest  event  is  underestimated  by  about  a 
factor  of  five. 

In  section  4,  we  have  converted  471  seismograms  from  103  earthquakes  in 
the  vicinity  of  Mammoth  Lakes,  California  into  equivalent  Wood-Anderson  seismo¬ 
grams  and  estimated  local  magnitude  Mi  following  the  original  definition  made 
by  Richter.  We  found  that  the  distance  correction  terms  given  by  Richter’s 
attenuation  curve  yield  systematically  larger  Mi  estimates  for  distances  less 
than  20  km.  ML  determined  at  near  distances  (  A<  10km  )  using  Richter’s  logA0 
values  are  as  much  as  1  magnitude  unit  greater  than  Mi  determined  at  regional 
distances  (  A >50  km  )  for  the  same  event.  This  result  is  similar  to  those  of 
recent  studies  utilizing  southern  California  and  northern  Baja  California  earth¬ 
quakes.  Seismograms  recorded  at  greater  distances  in  the  Great  Basin  (100  km 
<A<  600  km)  show  more  attenuation  (about  0.3  magnitude  units)  than  that 
predicted  by  Richter’s  curve.  We  have  determined  corrections  to  the  attenuation 
curve  which  removes  the  distance  dependence  observed  for  our  data.  We  also 
find  a  near  source  local  magnitude  saturation  for  events  with  Mi  S5.0  similar  to 
earlier  studies  with  California  earthquakes  and  explained  as  finite  source  size 
effects.  Using  the  revised  magnitude  scale  and  seismic  moments  estimated 
from  spectral  analysis  we  found  that  our  data  are  well  fit  by  the  straight  line 

log  M0  =(1.22±0.04  )ML  +  (16.98±0.15) 

for  1.5  UMi  ^6,3.  Comparison  of  our  data  to  predicted  M0  versus  Mi  values  sug¬ 
gest  a  constant  stress  drop  of  approximately  100  bars  for  events  with  Mi  <5.0 
and  increasingly  greater  stress  drop  for  the  larger  events. 


INFERENCES  ON  THE  CRUST  AND  UPPER  MANTLE 
STRUCTURE  OF  THE  CENTRAL  NEVADA  TEST  SITE  FROM 
TELESEISMIC  P-WAVE  RESIDUALS  AND  ITS 
IMPLICATIONS  FOR  THE  FAULTLESS  MAGNITUDE 
VERSUS  YIELD  ANOMALY 


David  E  Chavez  and  Keith  F  Priestley 


The  underground  nuclear  explosion  FAULTLESS  was  detonated  at  the  Cen¬ 
tral  Nevada  Test  Site  (CNTS)  in  Hot  Creek  Valley,  Nevada,  approximately  150  km 
north-northwest  of  the  Nevada  Test  Site  (NTS).  The  ISC  body  wave  magnitude, 
mb,  of  6.3  is  more  than  0.3  units  (a  factor  of  two  in  amplitude)  greater  than  that 
predicted  from  emperical  magnitude-yield  curves  for  NTS  explosions,  given  the 
announced  maximum  yield  of  1200  kilotons.  CNTS  is  physiographically  very 
similar  to  Pahute  Mesa  (NTS),  so  in  view  of  the  close  proximity  of  th  two  test 
sites  it  might  be  expected  that  the  magnitude-yield  curves  for  NTS  would  also 
apply  to  CNTS. 

Spence  (1974)  used  travel-time  residuals  for  NTS  explosions  recorded  at 
teleseismic  distances  to  infer  the  presence  of  a  high  velocity  structure  in  the 
mantle  beneath  the  14  m.y.  old  Silent  Canyon  caldera  at  Pahute  Mesa.  He  attri¬ 
buted  the  higher  velocities  to  material  depleted  of  its  partial  melt  component  by 
differentiation  and  by  eruption  through  the  caldera.  Later  studies  have 
confirmed  his  initial  observations.  In  order  to  see  if  any  anomalous  structures 
of  a  similar  nature  exist  beneath  CNTS  which  could  produce  the  discrepancy  in 
observed  mb  for  FAULTLESS,  we  have  examined  teleseismic  P-wave  residuals  as 
computed  from  a  seven-station  temporary  network  in  Hot  Creek  Valley. 

Figure  1  is  a  map  giving  the  generalized  geology,  station  locations,  and  the 
FAULTLESS  shot  point.  As  seen  from  the  map,  the  area  around  Hot  Creek  Valley 
is  dominated  by  a  series  of  nested  calderas  which  range  in  age  from  25  to  35 
m.y.;  these  calderas  are  10  to  20  m.y.  older  than  Silent  Canyon  caldera  at 
Pahute  Mesa.  However,  Hot  Creek  Valley  is  also  the  site  of  very  recent  volcan- 
ism.  For  instance,  much  of  the  Quaternary  basalt  shown  on  the  map  is 
sufficiently  young  that  no  erosion  is  evident. 

We  first  searched  for  local  lateral  variations  in  structure  by  computing 
teleseismic  P-wave  residuals  relative  to  the  average  residual  across  the  network 
for  each  event.  In  Figure  2  we  present  maps  which  show  relative  residual  as  a 
function  of  azimuth  for  each  CNTS  station.  Each  station  location  is  indicated  by 
a  circle  whose  radius  is  0.5  second.  The  circumference  represents  a  zero  resi¬ 
dual,  and  lines  drawn  from  the  circle  are  scaled  according  to  the  size  of  the  resi¬ 
dual.  Lines  inside  the  circle  are  negative  residuals  (early  arrivals)  and  lines  out¬ 
side  the  circle  are  positive  residuals  (late  arrivals).  The  azimuth  of  each  line 
corresponds  to  the  azimuth  of  approach  of  the  wavefront. 
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Beneath  each  map  is  an  array  diagram,  which  is  a  plot  of  mislocation  vec¬ 
tors,  i.e.,  the  vector  difference  between  the  expected  location  vector  (azimuth  of 
approach  and  slowness)  and  that  which  is  measured  by  fitting  a  plane  wave  to 
the  data.  The  existence  of  mislocation  vectors  is  due  to  differences  between  the 
earth  model  used  to  determine  the  expected  location  vector  and  the  real  earth. 
Mislocation  vectors  of  differing  magnitude  and/or  direction  are  indicative  of 
lateral  heterogeneities  between  the  source  and  receiver. 

The  first  pair  of  figures  are  for  uncorrected  data.  The  large  positive  residu¬ 
als  at  CNBC  are  due  to  approximately  1  km  of  sediments.  In  the  second  pair,  the 
data  have  been  corrected  for  surface  geology  which  results  in  a  systematic 
mislocation  to  the  southwest  as  evident  from  the  array  diagram.  This  misloca¬ 
tion  can  be  removed  by  correcting  the  data  for  a  Moho  striking  //  43°  £  with  3° 
dip.  This  interpretation  is  consistent  with  observations  of  Pn  velocity  across  the 
network  and  with  the  regional  Bouguer  gravity  anomaly. 

In  Figure  3,  we  give  x-y  plots  of  residual  versus  azimuth  for  the  data  which 
have  been  corrected  for  the  sediments  and  inferred  Moho  dip.  As  seen,  much  of 
each  observed  residual  is  removed  by  these  corrections.  The  late  arrivals  at 
CNBC  are  probably  due  to  insufficient  sediment  correction.  Two  stations,  CNPS 
and  CNEC,  continue  to  display  an  anomalous  pattern  in  that  for  events  arriving 
from  the  northwest  they  both  show  a  large  variation  in  residual.  For  emphasis, 
we  have  isolated  the  plots  for  these  two  stations  in  a  box  in  the  lower  right-hand 
corner  of  the  figure.  As  shown  by  the  two  plots  in  the  extreme  lower  right,  the 
residuals  at  CNPS  and  CNEC  show  a  very  definite  dependence  on  angle  of 
incidence.  The  anomaly  at  CNPS,  which  is  also  somewhat  evident  at  CNHC,  has 
not  yet  been  interpreted.  We  believe  that  the  anomaly  at  CNEC  is  related  to  its 
proximity  to  the  Lunar  Craters  volcanic  field. 

We  have  used  a  simple  geometrical  argument  to  qualitativly  estimate  the 
position  of  the  suspected  body  responsible  for  the  CNEC  anomaly.  Figure  4  is  a 
vie  r  along  a  line  passing  through  stations  CNHM  and  CNEC.  The  stippled  regions 
are  what  is  sampled  by  our  teleseismic  data,  where  the  heavy  shading  indicates 
normal  arrivals,  and  the  lighter  indicates  the  largest  possible  area  causing  the 
delay  at  CNEC.  We  observe  arrivals  as  late  as  0.3  sec,  comprable  to  that 
reported  by  Steeples  and  Iyer  (1976)  for  stations  in  Long  Valley  caldera  which 
they  attributed  to  anomalously  hot  rock.  Subsequent  studies  in  Long  Valley 
have  revealed  areas  of  S-wave  shadowing  associated  with  the  delay  that  have 
been  interpreted  as  evidence  for  the  presence  of  magma  (Sanders,  1984).  In 
view  of  the  young  volcanics  present  in  Lunar  Craters  it  is  not  unreasonable  to 
conclude  that  a  similar  situation  exits  here. 

Residuals  relative  to  event  average,  such  as  we  have  done,  will  not  reveal 
any  regional  velocity  anomalies,  so  we  have  also  computed  residuals  relative  to  a 
regional  event  average  using  many  stations  throughout  the  Great  Basin.  The 
results  are  shown  in  Figure  5a.  We  see  that  Hot  Creek  Valley  is  relatively  fast, 
suggesting  that,  as  in  Silent  Canyon,  a  high  velocity  vestige  of  a  root  to  the  Hot 
Creek  calderas  exits  at  depth.  Comparing  travel-time  residuals  the  explosion  in 
Hot  Creek  Valley  and  Silent  Canyon  relative  to  an  explosion  at  Yucca  Flat  (Figure 
5b)  indicates  slower  relative  velocities  beneath  CNTS  than  at  NTS.  .An  explana¬ 
tion  for  this  is  that  the  body  under  CNTS  is  probably  older  and  thus  has  had 
more  time  to  be  assimilated  into  the  surrounding  mantle. 

In  view  of  these  observations,  we  felt  that  ray  tracing  may  shed  light  on  the 
m.),  versus  yield  discrepancy  between  Hot  Creek  Valley  and  Pahute  Mesa.  Shown 
in  Figure  6  are  ray  paths  through  the  postulated  high  speed  bodies  beneath  both 
regions.  The  diagrams  indicate  that  for  Pahute  Mesa  events  there  is  a  caustic 
near  50  degrees  and  defocussing  at  greater  distances.  The  model  for  Hot  Creek 
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Valley  has  only  a  limited  effect. 

Most  mb  measurements  of  Nevada  nuclear  explosions  are  made  from  Euro¬ 
pean  stations  which  lie  in  the  Pahute  Mesa  shadow  zone.  Consequently,  we 
suspect  that  mb  vs  yield  averages  based  on  Pahute  Mesa  events  may  be  biased 
and  thus  when  applied  to  FAULTLESS  give  a  too  high  yield.  For  the  remainder  of 
the  contract  period  we  will  examine  WWSSN  recordings  of  FAULTLESS,  DUMONT, 
and  GREELEY  to  see  if  any  such  shadowing  effects  can  be  observed. 
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Figure  captions 


Figure  1.  Generalized  geologic  map  of  the  Hot  Creek  Valley,  Nevada  area  includ¬ 
ing  station  locations  and  the  FAULTLESS  shot  point.  Geology  from  Stewart 
and  Carison  (1977),  caldera  boundaries  from  Stewart  (1980). 


Figure  2.  Maps  of  travel  time  residuals  versus  azimuth  and  array  diagrams  (for¬ 
mat  described  in  the  text). 


Figure  3.  Plots  of  resiudal  versus  azimuth  for  all  stations  after  corrections  for 
surface  geology  and  inferred  Moho  dip  have  been  applied.  The  data  out¬ 
lined  in  the  lower  right  corner  are  for  stations  CNPS  and  CNEC  and  include 
plots  of  residual  versus  angle  of  incidence. 


Figure  4.  Sectional  view  along  a  line  passing  through  stations  CNHM  and  CNEC. 
Shading  described  in  the  text 


Figure  5.  a)  Map  of  residual  relative  to  event  average  which  includes  data  from 
several  Great  Basin  seismic  stations,  as  well  as  those  shown  on  the  map.  b') 
Map  of  residuals  for  FAULTLESS  and  GREELEY  relative  to  DUMONT  using  ISC 
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arrival  times. 


Figure  6.  a)  Ray  paths  through  the  high  speed  body  under  Pahute  Mesa  as  postu¬ 
lated  by  Spence  (1974).  b)  Ray  paths  through  the  high  speed  body  under 
Hot  Creek  Valley  as  suggested  by  our  data.  Figure  6.  a) 
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ABSTRACT 

Spectral  amplitudes  for  24  events  of  the  Mammoth  Lakes 
earthquake  sequence  have  been  determined  for  the  frequency 
range  0.1  to  10.0  Hz.  The  Mi  6  earthquake  that  occurred  at  1450 
UT  on  May  27,  1980  is  of  particular  interest  because  of  the  contro¬ 
versy  surrounding  its  mechanism.  A  comparison  of  spectral  levels 
determined  from  analysis  of  locally  recorded  strong-motion  and 
broad-band  digital  data,  and  spectral  levels  from  regionally 
recorded  surface  waves  extrapolated  back  to  the  source  yield  con¬ 
sistent  results  indicating  flat  spectra  in  the  band  0.1  -  1.0  Hz.  The 
specta  observed  in  this  study  do  not  show  pronounced  spectral 
peaks  predicted  by  theoretical  studies  of  gound  motion  due  to  the 
jerky  extension  of  a  fluid  driven  tensile  crack  containing  a  low 
viscosity  fluid.  We  have  found  nothing  in  the  spectra  of  the  Mi  6 
event  that  occurred  at  1450  UT,  on  May  27,  1980  nor  in  the  spectra 
of  the  aftershocks  to  distinguish  them  from  spectra  of  "tectonic" 
earthquakes.  However,  the  spectra  themselves  do  not  distinguish 
between  various  possible  explanations  for  the  non-double  couple 
source  mechanism  observed  in  moment  tensor  inversion  and  first 
motion  data  for  the  largest  events  of  the  Mammoth  Lakes  earth¬ 
quake  sequence 


Introduction 

The  Mammoth  Lakes  earthquakes  are  located  at  the  intersection  of  the 
Sierra  Nevada  frontal  fault  system  and  the  Long  Valley  caldera.  The  caldera  was 
formed  0.7  million  years  ago  by  collapse  and  subsidence  associated  with  the 
eruption  of  the  Bishop  tuff  (Bailey  ef  at.,  1976),  and  volcanism  has  continued  on 
a  reduced  scale  into  the  Holocene.  The  dominant  fault  in  the  area  south  of  the 
caldera  is  the  Hilton  Creek  fault  along  which  several  hundred  meters  of  pure 
normal  fault  displacement  has  occurred  since  the  formation  of  the  caldera  (Bai¬ 
ley  et  al.,  1976).  Extension  of  the  Hilton  Creek  fault  into  the  caldera  appears  to 
have  occurred  as  recently  as  0.3  million  years  ago,  possibly  indicating  that  the 
caldera  had  cooled  sufficiently  by  that  time  to  support  stresses  large  enough  to 
generate  earthquakes. 
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The  Mammoth  Lakes  earthquake  sequence  clearly  began  with  an  Ml  5  8 
event  on  October  4,  1978,  located  beneath  the  Wheeler  Crest  30  km  northwest  of 
Bishop,  California.  Following  this  event  activity  gradually  spread  to  the  north 
and  west.  The  main  energy  release,  which  began  on  May  25,  1980,  was  preceded 
for  over  a  year  by  Ml  3.5— 4.5  earthquake  swarms  in  the  region  of  the  Ml  6 
events.  These  swarms  generally  clustered  near  the  southern  boundary  of  the 
caldera,  and  migrated  eastward  towards  the  eventual  epicenter  of  the  first  Ml  6 
event  (Cramer  and  Toppozada,  1980).  The  main  episode  of  strain  release  com¬ 
menced  with  an  Ml  6. 1  (all  local  magnitude  estimates  for  the  main  events  are 
from  Uhrhammer  and  Ferguson,  1980)  event  at  1633  UT  on  May  25  (  hereafter 
referred  to  as  5/25/80--1633  ),  and  was  followed  by  a  second  Ml  6.0  event  at 
1649  UT.  These  initial  events  occurred  along  an  east-west  trend  near  the  south¬ 
ern  boundary  of  the  caldera.  Seismicity  then  migrated  southward  from  the  cal¬ 
dera  boundary.  Cramer  and  Toppozada  (1980)  note  that  no  events  at  the  magni¬ 
tude  5  level  occurred  south  of  the  southern  boundary  of  the  caldera  until  after 
the  Ml  6.1  event  at  1944  UT,  which  occurred  six  to  seven  kilometers  south  of  the 
first  Ml  6  event.  The  last  of  the  large  events  (Ml  6.2)  occurred  about  nine 
kilometers  south  of  the  1944  event,  at  1450  UT  on  May  27  (  hereafter  referred  to 
as  5/27/80—1450  ).  Lide  and  Ryall  (1984)  have  found  that  in  the  two  months  fol¬ 
lowing  the  main  events,  aftershocks  and  continuing  swarms  occurred  along  the 
southern  boundary  of  the  caldera,  and  along  several  NNE-SSW  trends  nearly 
orthogonal  to  the  caldera  boundary.  The  first  and  second  Ml  6  events  occurred 
along  the  first  of  these  trends,  and  the  first,  third  and  last  of  the  ML  6  events 
occurred  along  the  latter  trend.  The  epicenters  of  nearly  all  of  the  aftershocks 
and  all  of  the  Ml  6  events  are  well  to  the  west  of  the  projection  of  the  Hilton 
Creek  fault  at  depth  indicating  that  this  is  not  the  causative  fault  for  the  earth¬ 
quake  sequence. 

Considerable  controversy  has  arisen  regarding  the  source  mechanisms  of 
the  larger  events  of  the  Mammoth  Lakes  sequence.  Focal  mechanisms  derived 
from  Local  and  regional  first  motion  data  indicate  strike-slip  motion  with  a  NE- 
SW  T-axis  and  NW-SE  P-axis  (Cramer  and  Toppozada,  1960;  Ryall  and  Ryall,  1980). 
This  implies  right-lateral  motion  on  east-west  trending  faults,  or  left-lateral 
motion  on  north-south  trending  faults.  Unconstrained  moment-tensor  inversion 
of  long-period  teleseismic  body  and  surface-wave  data  result  in  a  large,  non- 
double-ccuple  component  for  the  first  and  last  of  the  Ml  6  events,  and  the  inver¬ 
sion  of  the  teleseismic  waveform  data  plus  teleseismic  first  motion  data  indicate 
mechanisms  for  the  first,  third  and  fourth  of  the  Ml  6  events  which  differ 
significantly  from  strike-slip  (Given  et  ai,  1982;  Barker  and  Langston,  1983). 

The  observed  non-double  couple  mechanism  may  arise  in  a  number  of  ways. 
Given  et  al.  (1982)  have  constrained  the  moment  tensor  to  be  a  double  -  couple 
and  interpreted  the  results  in  terms  of  left-lateral  oblique  slip  on  NNE  trending 
planes  dipping  eastward.  They  suggest  that  the  discrepancy  between  the 
regional  and  teleseismic  mechanisms  could  arise  from  either  structural  effects 
which  distort  the  teleseismically  observed  radiation  pattern,  or  complex  rupture 
during  the  event.  Distortion  of  the  outgoing  rays  may  result  from  complex 
structure  beneath  the  Long  Valley  caldera,  for  example  the  presence  of  magma 
chambers  which  were  first  proposed  by  Hill  (1976).  Similar  anomalies  in  focal 
mechanisms  have  been  reported  for  mid-oceanic  ridge  earthquakes,  and  these 
have  been  attributed  to  defocusing  of  the  ray  paths  due  to  magma  beneath  the 
rise  crest  (Solomon  and  Julian,  1974).  Wallace  (1984)  and  Lide  and  Ryall  (1984) 
have  interpreted  the  non-double  couple  mechanism  as  arising  from  a  multiple 
source.  Reversal  of  polarity  observed  on  short-period  and  long-period  instru¬ 
ments  at  the  same  site  (Wallace  et  at.,  1982)  suggest  a  complex  source  time 
function  with  the  short  period  arrivals  representing  the  failure  of  an  asperity 


-  3  - 


(possibly  with  a  strike  slip  mechanism,  and  the  long-period  arrivals  due  to  the 
overall  response  of  the  faulting  episode.  Similar  non-double  couple  mechanisms 
for  events  in  Greece,  have  been  shown  to  be  the  result  of  events  with  rupture  on 
multiple  planes  (D.  MacKenzie,  personal  communications,  1984). 

On  the  other  hand,  Julian  (1983)  has  suggested  that  the  results  of  inverting 
the  surface  wave  and  long-period  body  wave  data,  and  taking  the  regional  and 
teleseismic  first  motion  data  as  a  whole,  are  better  fit  by  a  compensated  linear 
vector  dipole  (CLVD)  mechanism  (KnopofT  and  Randall,  1970),  than  they  are  by  a 
double-couple  mechanism.  Julian  has  noted  that  the  CLVD  mechansm  is  con¬ 
sistent  with  all  but  four  and  two  of  the  first  motion  data,  respectively,  for  the 
5/25/80—1633  and  5/27/80—1450  Ml  6  events,  as  opposed  to  32  and  16  incon¬ 
sistent  data  for  the  double  couple  mechanism,  and  has  suggested  that  the 
5/25/80—1633  and  5/27/80—1450  events  are  due  to  injection  of  a  fluid  into  an 
expanding  crack,  i.e.,  dike  formation.  Aki  (1984),  utilizing  the  calculations  of 
Chouet  (1981),  has  suggested  that  the  non-double  couple  mechanism  results 
from  the  rapid  opening  of  a  channel  between  two  pre-existing  cracks  filled  with 
magma. 

In  this  study  we  compare  spectral  amplitudes  of  24  Mammoth  Lakes  earth¬ 
quakes  ranging  in  magnitude  from  Ml  3.5  to  Mi  6.2,  over  the  frequency  band  ss 
0.05  -  20  Hz.  Spectral  amplitudes  from  approximately  0.5  to  20  Hz  are  derived 
from  spectral  analysis  of  locally  recorded  body-wave  data.  Spectral  amplitudes 
from  approximately  0.05  to  0.125  Hz  are  derived  from  analysis  of  regionally 
recorded  surface  wave  data.  The  wide  bandwidth  of  the  spectral  amplitude 
measurements  allows  us  to  compare  the  observed  spectral  shapes  with  theoreti¬ 
cal  spectral  shapes  for  ground  motion  resulting  from  "tectonic”  earthquakes 
and  ground  motion  predicted  for  a  propagating  tensile  crack. 

Data 

This  study  is  based  on  seismograms  of  the  Mammoth  Lakes  earthquakes 
recorded  on  the  seismic  stations  shown  in  Figure  1.  The  events,  shown  as  dots, 
comprise  only  a  small  fraction  of  the  total  number  of  events  which  have 
occurred  in  the  Mammoth  Lakes  area  since  October,  1978.  Locations  after  May 
27,  1980  are  largely  from  Lide  (written  communication,  1984).  Locations  for 
some  of  the  earlier  events  are  from  the  University  of  Nevada  or  University  of 
California  at  Berkeley  bulletins.  The  events  selected  for  analysis  were  those 
recorded  on  at  least  one  local  strong-motion  or  digital  seismograph  and  one 
long-period  regional  seismograph.  The  SMA-1  strong-motion  accelerographs, 
shown  as  triangles  in  Figure  1,  were  installed  by  the  California  Division  of  Mines 
and  Geology  (Turpin,  1980),  or  by  the  University  of  Southern  California  (USC) 
(Moslem  et  a l,  1983).  We  also  used  shear  wave  spectra  supplied  by  Ralph 
Archuleta  from  data  recorded  by  the  U.S.  Geological  Survey.  Digital  stations, 
shown  by  squares  in  Figure  1  and  described  by  Archuleta  et  at.  (1982),  consist  of 
either  force-balance  accelerometer  or  velocity  transducers  and  Sprengnether 
DR-100  recorders.  For  events  along  the  southern  boundary  of  the  caldera,  the 
SMA-1  and  digital  seismographs  cover  a  large  fraction  of  the  focal  sphere  com¬ 
pared  with  many  previous  studies  of  a  similar  nature.  The  radiation  pattern  of 
events  south  of  the  caldera  are  not  as  well  sampled. 

The  inset  map  in  Figure  1  shows  the  locations  of  the  regional  long-period 
stations  from  which  seismograms  were  used  to  determine  surface  wave  magni¬ 
tude.  These  stations  completely  surround  the  Mammoth  Lakes  region  and,  are 
all  WWSSN  stations,  with  the  exception  of  the  California  Institute  of  Technology 
(CIT)  station  at  Pasadena  (PAS). 
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Analysis 

Determination  of  Body-Wave  Moment:  Displacement  amplitude  spectra  were 
computed  by  Fourier  transforming  the  S-wave  arrival  on  the  strong-motion 
accelerograms,  and  then  dividing  the  acceleration  spectra  by  cj“2  to  give  the  dis¬ 
placement  spectra.  The  time  window  was  chosen  to  be  long  enough  to  include 
the  entire  S-wave  arrival,  but  short  enough  to  minimize  the  inclusion  of  surface- 
wave  energy.  The  spectra  have  been  normalized  to  a  distance  of  10  km.  Since 
the  precise  locations  of  the  foci  of  the  events  are  not  known,  we  have  not  used 
the  epicenter  and  seismograph  coordinates  to  determine  the  distance,  but  have 
estimated  the  epicentral  distance  R  from  the  S-P  times  assuming  an  avarage  S- 
wave  velocity  of  3.4  km  s_l  and  an  average  P-wave  velocity  of  5.8  km  s~{.  In 
addition,  we  have  not  rotated  the  data  into  transverse  and  radial  components, 
but  have  computed  the  spectra  for  both  horizontal  components  and  computed 
the  vector  sum  of  these  two  values. 


The  spectra  were  interpreted  in  terms  of  a  flat,  low-frequency  asymptote 
n0.  a  high  frequency  falloff.  and  the  corner  frequency  fc  at  which  the  two 
trends  intersect.  The  far-field  spectral  parameters  were  then  related  to  the 
source  parameter  seismic  moment  M0  using  the  relationship  given  by  Brune 
(1970,  1971) 


_  4rrp/g3noiP 
M°~  kR ^ 


(1) 


where  p  is  the  density  (p  =  2.8  g  cm-3),  @  is  the  shear-wave  velocity 

(jS  =  3.5  km  s-1),  R$v  is  the  RMS  average  of  the  radiation  pattern  (taken  here  as 
0.6),  and  k  is  a  correction  factor  for  the  amplification  of  SH  waves  at  the  free 
surface  (taken  as  2).  The  correction  of  the  amplitude  of  the  S-wave  on  being 
reflected  at  the  free  surface  assumes  that  the  waves  are  horizontally  polarized 
{SH).  Since  we  have  not  rotated  the  data  into  transverse  and  radial  com¬ 
ponents,  the  data  are  in  general  some  mixture  of  SH  and  ST  motion.  Abe 
(1974)  discusses  experimental  studies  of  amplification  on  reflection  and  con¬ 
cludes  that  using  k  =  2  applies  approximately  to  SH-wave  signals  containing 
some  ST  motion.  Helmberger  and  Malone  (1975)  discuss  problems  due  to  multi¬ 
ple  arrivals  caused  by  reflections  in  deducing  source  parameters  from  S-waves 
recorded  at  epicentral  distances  much  greater  than  the  source  depth.  The 
station-source  configuration  shown  in  Figure  1  indicates  that  this  should  not  be 
a  problem  in  this  study. 

A  source  of  error  inherent  in  the  estimate  of  n0  and  consequently  in  the 
estimates  of  Ma  comes  from  the  assumption  of  a  single,  average  correction  for 
radiation  pattern.  If  sufficient  instruments  were  available,  or  if  the  source 
parameters  including  direction  and  velocity  of  rupture  were  known,  this  source 
of  error  could  be  reduced.  A  value  of  0.6  was  used  for  R$9  (Thatcher  and  Hanks, 
1973).  At  some  azimuths,  the  amplitude  could  be  1.67  times  larger  than  the 
amplitude  determined  assuming  this  correction,  while  near  a  node  in  the  radia¬ 
tion  pattern  the  amplitude  could  be  quite  small.  Furthermore,  if  the  source 
dimension  is  comparable  to  the  wavelength  involved,  source  propagation  might 
distort  the  radiation  pattern  by  focusing  or  defocusing  energy,  thereby  changing 
the  apparent  corner  frequency  for  different  directions  and  possibly  reducing  the 
calculated  moment. 

No  Q  correction  was  made  in  determining  the  spectra.  Archuleta  et  al. 
(1982)  found  that  whole-path  Q  had  a  small  effect  on  the  measurement  of  H0  for 
Mammoth  Lakes  events. 
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Determination  of  Surface  Wave  Moment:  Tucker  and  Brune  (1977)  determined 
the  seismic  moment  of  seven  aftershocks  of  the  San  Fernando  earthquake  from 
the  surface  wave  spectra.  This  gave  an  estimate  of  seismic  moment  indepen¬ 
dent  of  that  obtained  from  S-wave  spectra  and  one  based  on  lower  frequencies 
(0.03  -  0.10  Hz  in  their  case).  They  made  a  careful  estimate  of  the  errors 
involved  in  estimating  the  seismic  moment  from  the  surface  wave  spectra.  The 
errors  included  procedural  errors,  path  structure  and  attenuation,  lateral 
refraction  (multipathing),  source  finiteness,  and  uncertainties  in  fault  parame¬ 
ters  (depth,  dip,  strike,  and  slip  vector)  in  correcting  for  radiation  pattern. 
They  concluded  that  these  effects  could  lead  to  a  maximum  error  in  the  final 
result  of  no  more  than  a  factor  of  4. 

In  this  study  we  also  determined  the  seismic  moment  from  surface  waves 
for  comparison  with  seismic  moments  determined  from  the  S-wave  spectra. 
However,  we  used  a  different  procedure  which  allowed  us  to  determine  the 
seismic  moment  from  a  large  number  of  events  without  the  time-consuming 
digitization  of  the  surface  wave  seismograms. 

Surface  wave  magnitude  (Ms)  defined  by  the  formula  (Marshall  and  Basham. 
1972) 

Ms  =  Log  A  +  fl'(A)  +  P(T)  (2) 

was  determined  from  seismograms  of  long-period  seismograph  stations,  pri¬ 
marily  in  the  western  United  States.  In  equation  (2),  A  is  one-half  the  peak-to- 
peak  ground  amplitude  of  the  maximum  amplitude  Rayleigh  wave  with  period  T 
at  distance  A.  B'( A)  corrects  for  average  attenuation,  scattering,  geometrical 
spreading  and  refraction;  and  P(T)  is  a  correction  factor  for  dispersion  and 
adlows  measurement  of  the  surface-wave  amplitude  at  any  period.  Marshall  and 
Basham  (1972)  have  tabulated  P(T)  for  continental  North  America  for  periods 
ranging  from  10  to  40  seconds.  In  some  cases  the  maximum  surface-wave  ampli¬ 
tude  observed  for  the  Mammoth  Lakes  earthquakes  was  in  the  8  to  10  second 
range.  Since  P(T)  is  a  smoothly  varying  function,  we  have  extrapolated  P(T)  to 
eight  seconds  period.  The  surface  wave-magnitudes  for  the  Mammoth  Lakes 
events  studied  here  are  tabulated  in  Table  1. 

The  surface-wave  moment  is  linearly  related  to  the  surface  wave  amplitude, 
and  thus  to  the  surface  wave  magnitude.  We  have  determined  a  relationship 
between  the  surface-wave  magnitude  and  the  seismic  moment  in  the  following 
manner.  Surface  wave  synthetics  were  generated  by  a  computer  routine  (writ¬ 
ten  by  Guy  Masters)  for  an  earthquake  with  M0  =  10zl  dyne-cm  at  a  number  of 
regional  seismographs  from  which  data  were  used  to  determine  the  surface-wave 
magnitudes.  From  the  synthetic  surface  wave  seismograms,  we  determined  the 
surface  wave  magnitude  in  the  same  manner  as  discussed  above  for  the  Mam¬ 
moth  Lakes  earthquakes,  and  thus  determined  a  relationship  between  the  sur¬ 
face  wave  magnitude  and  seismic  moment. 

Factors  affecting  the  calculation  of  the  synthetic  surface  wave  seismograms 
are  the  earth  velocity  and  attenuation  structure  and  the  fault  parameters  of  the 
individual  earthquakes.  Since  the  events  are  located  at  the  western  edge  of  the 
Great  Basin,  while  most  of  the  seismograph  stations  are  located  in  western  North 
America,  we  have  used  the  velocity  structure  determined  for  the  Great  Basin 
from  surface  wave  dispersion  measurements  (Priestley  and  Brune,  1978).  The 
attenuation  structure  for  the  Great  Basin  by  Patton  and  Taylor  (1984)  was  used, 
with  a  slight  lowering  of  Q0  from  that  of  Patton  and  Taylor  in  the  surficial  layer. 
Since  we  do  not  have  fault-plane  solutions  nor  do  we  know  the  focal  depth  for  the 
majority  of  events  we  have  studied,  we  determined  an  average  relationship  for 
the  Mammoth  Lakes  area  in  the  following  manner.  Vetter  and  Ryall  (1983)  have 
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determined  fault  plane  solutions  for  a  large  number  of  events  in  the  Mammoth 
Lakes  area,  and  found  that  the  mechanisms  can  generally  be  divided  into  two 
groups.  The  first  group  consists  of  earthquakes  with  strike-slip  mechanisms  ( 
average  plane  of  strike  N14°E,  dip  80°ESE,  slip  N16°E  or  strike  N73°W,  dip 
80°SSW,  slip  N76°E  )  which  are  most  commonly  at  depths  less  than  9  km.  The 
second  group  consists  of  earthquakes  with  oblique  or  normal  mechanisms  ( 
average  planes  of  strike  N03°E,  dip  64°E,  slip  N40°E  or  strike  N5Q°W,  dip  38°SW, 
and  slip  N86°W  )usuallv  at  depths  greater  than  9  km.  Synthetic  surface  wave 
seismograms  were  calculated  for  both  mechanisms  for  2.  5,  8,  11  and  14  km 
depth.  Surface-wave  magnitudes  were  determined  from  the  10  records  at  each 
station  and  an  average  conversion  factor  determined  for  both  focal  mechanisms 
and  all  focal  depths.  The  resulting  relationship  between  surface  wave  magnitude 
Ms  and  seismic  moment  Mo  for  the  Mammoth  Lakes  area  is  given  by  the  equa¬ 
tion 

log  M0  =  Ms  +  19.40(x0.06)  .  (3) 

To  verify  the  validity  of  our  moment-surface  wave  magnitude  relationship, 
we  can  compare  equation  (3)  with  other  moment  magnitude  relations  esta¬ 
blished  for  the  western  United  States.  Wyss  and  Brune  (1968)  originally  esta¬ 
blished  a  moment-local  magnitude  relationship  for  events  less  than  magnitude  8 
in  the  western  U.S.  Seismic  moments  were  determined  for  272  events  which 
averaged  over  a  wide  range  of  tectonic  regions  (  California,  Nevada,  Arizona, 
Utah  and  Baja  California  ).  These  were  found  to  be  best  fit  by  the  curve  given  by 
the  equation 

log  Mo  =  1.7  Ml  +  15. 1  3  <  Ml  <  6 

Using  the  Ms  -  ML  relationship  given  by  Wyss  and  Brune  (1968),  this  becomes  in 
terms  of  Ms(G )  (the  surface  wave  magnitude  defined  by  Gutenberg  (1945)) 

log  M0  =  MS(C)  +  19.2  3  <  Ms  <  6  . 

This  definition  can  be  compared  with  the  relationship  derived  directly  from 
Gutenoerg's  (1945)  definition  of  surface  wave  magnitude.  According  to  that 
definition,  a  magnitude  6  earthquake  produces  a  far-field  displacement  of 
100 /im  at  distance  22°  for  surface  waves  of  period  20  seconds.  In  terms  of 
moment  Gutenburg’s  relationship  becomes 

log  M0  =  MS(C)  +  19.3  . 

This  point  thus  represents  average  of  the  numerous  observations  on  which  the 
surface  wave  magnitude  was  based.  As  pointed  out  by  Richter  (1958),  the  scale 
was  adjusted  to  agree  with  the  local  magnitude  values  of  6  to  7.  Tucker  and 
Brune  (1977)  did  not  explicitly  give  a  relationship  between  surface  wave  magni¬ 
tude  and  seismic  moment.  However,  from  the  values  given  in  their  Tables  1  and 
3,  the  average  relationship  is 

log  M0  =  Ms  +  19.3  . 

This  relationship,  determined  using  a  Gutenburg  continental  earth  model  and 
the  Mars  hall -Basham  definition  of  surface  wave  magnitude,  is  in  good  agreement 
with  the  relationship  based  on  the  Gutenburg  definition  of  surface  wave  magni¬ 
tude.  Comparing  these  various  relationships  between  Ms  and  M0  we  conclude 
that  the  use  of  the  Great  Basin  velocity  (Priestley  and  Brune,  1978)  and  attenua¬ 
tion  (Patton  and  Taylor,  1984)  model  has  increased  the  estimated  moment  of  the 
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events  by  26%  over  that  determined  for  the  Gutenburg  continental  model. 

Equation  (3)  allows  us  to  rapidly  determine  the  seismic  moment  M q  from 
the  surface  wave  magnitude  observation  in  Table  1.  The  resulting  values  of 
seismic  moment  (see  Table  1)  were  converted  to  an  equivalent  value  of  Ho  by 
solving  equation  (1)  for  Go  and  substituting  the  seismic  moment  derived  from 
the  surface  wave  observations.  This  gives  values  of  G0  at  frequencies  of  the  sur¬ 
face  waves  for  comparison  with  the  spectra  obtained  from  the  body-wave  record 
on  the  accelerometers. 

Error  analysis:  Errors  in  fi0  (  and  consequently  in  M0  )  were  estimated  by  calcu¬ 
lating  an  average  spectral  value  and  a  standard  deviation  where  more  than  one 
recording  was  available  for  an  event.  The  average  values  of  fl0  was  determined 
by  the  formula 


1  NS 

i] 

<  G0  >  =  antxlog  ^  £  log 

Go^t/10] 

where  G0.  is  the  long-period  spectral  level  at  the  ith  station  ,  /?*  is  the  distance 
to  the  Ith  station,  and  NS  is  the  number  of  stations  recording  the  event.  The 
standard  deviation  of  the  log  average  value  of  no  was  determined  by  the  formula 


s.d.  (log  <  G0  >)  = 
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A  multiplicative  error  factor  is  then  determined  from 


EQ  o=  antilog 


s.d.  (log<G0  >) 


As  pointed  out  by  Archuleta  et  al.  (1902),  when  calculating  average  values  of  fl0 
or  its  derivative  M0,  it  is  necessary  to  compute  the  averages  using  the  equations 
in  the  above  forms  to  give  equal  weight  to  each  observation.  If  a  simple  arith¬ 
metical  average  value  is  determined  it  will  be  biased  toward  larger  values,  since 
the  errors  associated  with  Go  are  lognormally  distributed. 

Observed  Spectra:  Figures  2  and  3  are  displacement  spectra  for  two  events  of 
the  24  Mammoth  Lakes  earthquakes  we  studied.  All  spectra  are  computed  from 
strong-motion  accelerograms  of  the  USC-CIT  data  set.  The  plotted  spectra  are 
the  vectorial  average  of  the  two  orthogonal  components,  and  have  been  normal¬ 
ized  to  a  hypocentral  distance  of  10  km.  The  arrows  indicate  our  picks  of  the 
corner  frequency  and  long-period  spectral  level.  The  triangle  with  error  bars 
denotes  the  average  spectral  level  determined  from  the  surface  wave  measure¬ 
ments. 

Figures  2a-b  are  two  of  14  spectra  for  an  Mi  4.9  event  occurring  at  1516  UT 
on  May  31,  1980  (  hereafter  referred  to  as  5/31/80—1516  ).  The  focus  of  this 
event  is  nearly  identical  to  that  of  the  5/25/80-1633  event.  The  average  value 
of  the  long-period  spectral  level  <  G0  >  measured  from  the  14  spectra  is 
1.03x  10-1  cm-sec,  corresponding  to  M0-  1.34x  1023  dyne-cm.  Our  values  include 
esimates  of  G0  from  the  data  of  Archuleta  et  al.  (1982)  (six  observations),  plus 
eight  additional  spectra  from  the  USC-CIT  strong-motion  data.  Our  value  of 
<  G0  >  does  not  differ  significantly  from  their.  The  seismic  moment  from  the 
surface  wave  observation  is  1. 1  x  1023  dyne-cm  which  is  not  significantly  different 
from  the  body-wave  value. 

Spectra  of  the  largest  earthquake  we  studied,  the  5/27/80—1450  Mi  6 
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shock,  are  shown  in  Figure  3a-b.  This  event  is  of  particular  interest  because  of 
the  controversy  concerning  its  mechanism.  Four  values  of  moment  have  been 
published  for  this  event.  Archuleta  et  al.  (1982)  determined  a  moment  of 
2.33  x  1024  dyn  e-cm  from  the  spectral  analysis  of  two  locally  recorded  seismo¬ 
grams.  Uhrhammer  and  Ferguson  (1980)  estimated  a  moment  of  5.0xl024 
dyne-cm  from  the  spectra  of  51.2  seconds  of  the  broad-band  displacement 
record  of  the  vertical  component  seismogram  at  Jamestown,  California,  140  km 
west  of  Mammoth  Lakes,  Barker  and  Langston  (1983)  inverted  seven  long-period 
teleseismic  P-waveforms  and  five  long-period  teleseismic  SH-waveforms  to 
obtain  the  moment  tensor  for  this  event.  They  found  a  moment  of  1.03xl02S 
dyne-cm  and  36%  CLVD.  Given  et  al.  (1982)  inverted  long-period  surface  wave 
from  this  event  and  used  this  result  to  model  the  long-period  body-waves.  They 
found  a  moment  of  l.lOx  1025  dyne-cm.  Our  values  of  7.26  x  1024  and  7.20  x  IQ24 
dyne-cm  for  the  body-  and  surface-wave  observations,  respectively,  are  some¬ 
what  higher  than  the  values  reported  by  Archuleta  et  al.  (1982)  but  in  reason¬ 
able  agreement  with  the  other  reported  values  of  moment. 

Discussion 

Figure  4  is  a  plot  comparing  the  surface  wave  moments  and  body-wave 
moments  of  the  Mammoth  Lakes  earthquakes  from  Table  1.  The  reference  line 
through  the  points  has  unit  slope  and  intercept  zero.  Those  points  lying  above 
the  line  have  larger  body-wave  moments  compared  to  the  observed  surface  wave 
moment,  while  those  below  the  line  have  larger  surface-wave  moments  com¬ 
pared  to  the  observed  body-wave  moment.  Vertical  error  bars  denote  the 
spread  in  the  observed  body-wave  moment  while  horizontal  error  bars  denote 
the  spread  in  the  observed  surface-wave  moment.  It  is  clear  from  this  figure 
that  there  are  no  significant  differences  between  the  moments  determined  at  8 
to  20  seconds  from  surface  waves,  and  the  moments  determined  at  periods  near 
1  second  from  body-waves.  Where  we  have  analysed  events  in  common,  our 
results  from  the  surface  wave  measurements  support  the  moment  estimates  of 
Archuleta  et  al.  (1982)  and  thus  their  findings  of  nearly  constant  stress  drop  (Air 
~  50  bars)  for  events  with  seismic  moments  greater  than  approximately  1  x  1021 
dyne-cm.  All  of  the  events  we  have  studied  are  larger  than  1  x  1031  dyne-cm  and 
therefore  we  cannot  substantiate  the  dependence  of  stress  drop  on  seismic 
moment  which  Archuleta  et  al.  (1982)  observed  for  smaller  events.  However,  the 
wider  bandwidth  of  our  spectra  permit  us  to  address  the  problem  of  the  source 
mechanism  controversy  surrounding  the  Mammoth  Lakes  events. 

Aki  et  al.  (1977)  and  Chouet  (1981)  derive  theoretical  ground  motion  in  the 
far  and  near-feld  respectively,  due  to  the  extension  of  a  fluid-driven  tensile 
crack  embedded  in  a  layered  half-space.  In  both  studies,  the  source  was  a  jerky 
opening  of  a  channel  ahead  of  the  crack  tip  due  to  excess  pressure  of  the  fluid  in 
the  crack.  Seismic  waves  are  generated  by  the  vibrations  of  the  newly  created 
crack  walls  due  to  the  rapid  application  of  excess  fluid  pressure.  Both  the  far- 
field  and  near-field  studies  show  a  peaked  character  to  the  ground  motion  spec¬ 
tra.  in  contrast  to  "tectonic"  earthquake  spectra  which  are  flat  to  some  corner 
frequency  and  then  fall  off  at  higher  frequencies.  Chouet  (1981)  varied  the 
model  parameters  to  analyze  the  effects  of  fluid  compressibility,  source  depth 
and  structure  of  the  medium  on  the  ground  motion.  He  found  that  the  spectral 
peaks  depend  on  the  source  geometry,  the  medium  characteristics,  the  receiver 
position,  the  component  of  ground  motion  being  studied,  and  critically  on  the 

crack  stiffness  factor  C  defined  as  C  =  — - —  where  b  is  the  bulk  modules  of  the 

a  [A 

fluid,  L  is  the  full  length  of  the  crack,  d.  is  the  crack  width,  and  fj.  is  the  rigidity 
of  the  solid.  The  crack  vibrates  most  easily  when  the  crack  is  empty  (C  =  0)  (see 
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Aki  et  a!..,  1977,  Figures  6-9).  The  spectral  peak  becomes  broadened  however 
when  the  crack  stiffness  factor  is  about  5,  implying  that  the  bulk  modulus  of  the 
magma  is  low  enough  to  contain  bubbles  (Aki,  1984).  Chouet's  (1981)  model  cal¬ 
culations  show  a  pattern  similar  to  the  non-double  couple  mechanism  observed 
for  the  largest  Mammoth  Lakes  earthquakes  (5/25/80--1633,  5/27/80—1450). 
Utilizing  Chouet's  (1981)  calculations,  Aki  (1984)  has  proposed  that  these  events 
result  from  the  sudden  opening  of  a  channel  between  two  pre-existing  fluid-filled 
cracks,  due  to  a  higher  fluid  pressure  in  one  of  the  cracks.  The  pressure  in  the 
magma  drops  due  to  the  increased  total  crack  volume  and  the  crack  tip  closes 
as  a  result. 

However,  as  shown  in  Figure  4.  none  of  the  events  we  have  studied  show 
significant  differences  in  level  in  the  frequency  band  0.1  to  1.0  Hz  and  therefore 
do  not  show  the  spectral  peaking  predicted  by  the  Aki  (1984)  model  for  a  low 
viscosity  fluid.  Most  of  the  strong-motion  and  digital  recording  seismographs  for 
the  5/27/80—1450  event  were  located  to  the  north  of  this  event  thus  sampling 
only  a  portion  of  the  focal  sphere,  hence  it  may  be  argued  that  we  have  sampled 
that  part  of  the  radiation  pattern  where  the  peaked  nature  of  the  spectrum  is 
least  prominent.  Figure  5  shows  the  distribution  of  spectra  for  the  Mi  4.9 
5/31/80—1516  aftershock  which  occurred  in  almost  the  same  location  as  the 
5/25/80—1633  Mi  6  event  which  also  had  a  non-double  couple  mechanism  (Given 
et  at.,  1982;  Julian,  1983).  The  focal  mechanism  from  local  and  regional  first- 
motion  data  for  the  5/31/80-1516  event  is  from  Vetter  (written  communication, 
1984)  and  is  very  similar  to  the  solution  independently  determined  by  Archuleta 
et  al.  (1982).  Figure  5  shows  an  upper-hemisphere  plot  with  the  northwest  and 
southeast  being  the  compression  quadrants.  This  solution  is  similar  to  the  focal 
mechanism  given  by  Cramer  and  Toppozada  (1980)  and  Ryall  and  Ryall  (1980) 
from  Local  and  regional  first  motion  data  for  the  5/25/80—1633  event.  The  loca¬ 
tion  of  the  stations  from  which  the  spectra  were  computed  are  indicated  on  the 
focal  sphere.  Some  of  the  individual  spectra  show  a  suggestion  of  peaking  near  1 
Hz,  but  not  of  the  magnitude  predicted  by  Aki  et  at.  (1977)  and  Chouet  (1981). 
We  find  nothing  in  the  spectra  of  the  Mammoth  Lakes  events  we  have  studied 
which  distinguish  them  from  spectra  of  "tectonic"  events  observed  in  numerous 
previous  studies.  The  flat  nature  of  the  spectra  does  not  rule  out  the  opening  of 
a  tensile  crack  as  the  source  of  the  elastic  waves  for  these  earthquakes.  It  does 
preclude  the  opening  of  a  tensile  crack  with  the  subsequent  injection  of  a  fluid  of 
low  viscosity. 

Summary 

Spectra  for  24  Mammoth  Lakes  earthquakes  greater  than  1  x  1021  dyne-cm, 
including  the  Mi  6  event  occurring  at  1450  UT  on  May  27,  1980,  have  been  deter¬ 
mined.  The  short-period  portion  of  the  spectra  (0.5  -  10.0  Hz)  have  been  deter¬ 
mined  from  spectral  analysis  of  the  locally  recorded  strong-motion  and  broad¬ 
band  digital  data.  A  relationship  has  been  derived  between  the  Marshall-Basham 
surface-wave  magnitude  and  seismic  moment  allowing  us  to  determine  the 
seismic  moment  directly  from  time  domain  amplitude  measurements  of  the  sur¬ 
face  waves  recorded  at  regional  stations.  This  gives  an  estimate  of  the  seismic 
moment  independent  of  that  obtained  from  the  S-wave  spectra  and  one  based  on 
lower  frequencies  (0.05  -  0.125  Hz).  A  comparison  of  the  seismic  moments  indi¬ 
cate  that  for  any  one  event,  no  significant  difference  in  the  spectral  level  exist 
between  the  values  determined  at  surface  wave  frequencies  (0.05  to  0.125  Hz  ) 
and  the  S-wave  corner  frequency  (  1  to  5  Hz  ).  The  events  studied,  and  in  partic¬ 
ular  the  5/27/80-1450  ML  6  event,  do  not  show  peaked  spectra  predicted  from 
theoretical  analysis  of  the  ground  motion  resulting  from  a  fluid-driven,  pro¬ 
pagating  tensile  crack.  We  have  found  nothing  in  the  spectra  of  these  events  to 
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distinguish  them  from  "tectonic”  earthquakes.  The  absence  of  a  peaked  spec¬ 
trum  does  not  rule  out  dike  injection  as  the  source  of  these  earthquakes,  but  it 
does  rule  out  a  particular  model  of  magmatic  intrusion,  namely  that  of  Chouet 
(1981)  and  Aki  (1984),  for  a  low  viscosity  fluid.  As  Aki  (1984)  has  pointed  out,  the 
opening  of  a  channel  between  two  pre-existing  fluid  filled  cracks  would  Largely 
alleviated  this  problem  by  involving  the  movement  of  only  minor  amounts  of 
viscous  material,  which  could  be  accomplished  in  an  adequately  small  time. 
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Figure  Captions 

Figure  1.  Location  maps  showing  the  events  studied  and  the  stations  from  which 
seismograms  were  analysed.  In  the  lower  right-hand  map.  the  solid  circles 
denote  locations  of  events,  triangles  denote  location  of  strong-motion 
accelerographs,  and  squares  denote  location  of  digital  seismographs.  The 
upper  left-hand  man  shows  the  locations  of  the  long-period  seismographs 
used  for  the  surface  wave  analysis. 


Figure  2.  Spectra  from  strong-motion  recordings  at  the  sites  "McGee  Creek  Inn" 
and  "Long  Valley  Fire  Station"  for  the  Mi  4.9  event  occurring  at  1516  UT  on 
May  31,1980.  Spectra  are  the  vectorial  sum  of  the  two  orthogonal  com¬ 
ponents,  and  have  been  normalized  to  10  km.  The  long-period  spectral 
level  picked  from  these  spectra  are  denoted  by  the  arrow,  and  the  long- 
period  spectral  level  determined  from  the  surface-wave  observations  are 
denoted  by  the  triangle  with  error  bars. 


Figure  3.  Spectra  from  strong-motior}  recordings  at  the  sites  "Cash  Baugh 
Ranch"  and  "Long  Valley  Fire  Station"  for  the  Mi  6.2  event  occurring  at 
1450  UT  on  May  27, 1980.  The  notation  is  the  same  as  for  figure  2. 
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Figure  4.  Comparison  of  the  body-wave  and  surface-wave  determined  moments 
for  the  24  events  studied.  The  line  is  for  reference  and  is  of  unit  slope  and 
zero  intercept. 

Figure  5.  Variation  of  the  spectra  for  the  event  occurring  at  1516  UT  on  May  31, 
1980,  as  a  function  of  position  on  the  focal  sphere.  The  focal  mechanism  is 
from  Vetter  (written  communication,  1984). 
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Introduction 

Earthquakes  with  similar  short-period  amplitudes,  as  evidenced  by  similar 
local  magnitude  Mi,  can  show  large  differences  in  long-period  amplitudes.  These 
differences  can  arise  due  to  variations  in  the  properties  of  the  path,  the  source 
depth,  and  source  parameters.  Understanding  the  differences  in  the  relative 
excitation  of  high  frequency  (0.5-10.0  Hz)  and  low  frequency  (0.05-0.125  Hz) 
seismic  waves,  and  the  propagation  of  this  energy  at  regional  distance,  is  impor¬ 
tant  for  understanding  source  mechanisms,  and  more  specifically  for  evaluating 
the  possibilities  of  discriminating  small  explosions  from  earthquakes  at  regionai 
distance.  The  mb-Ms  discriminant,  which  compares  the  relative  excitation  of 
body  versus  surface  waves,  is  presently  one  of  the  primary  techniques  for  distin¬ 
guishing  explosions  from  earthquakes.  The  theoretical  basis  for  the  discrim¬ 
inant  is  still  not  well  understood,  it  is  thought  to  arise  from  differences  in  one  or 
several  of  the  following  source  parameters:  (l)  source  dimension,  (2)  source 
rupture  duration,  (3)  time  function  of  the  displacement  at  various  points  on  the 
fault  for  an  earthquake  or  at  the  elastic  radius  for  an  explosion,  and  (4)  source 
symmetry. 

Regional  variations  in  the  relative  excitation  of  high  and  low  frequency 
seismic  waves  generated  by  earthquakes  in  western  North  America  have  been 
noted  in  a  number  of  studies.  Brune,  Espinosa  and  Oliver  (1963)  found  a  wide 
range  in  the  surface  wave  excitation,  for  a  given  ML,  for  earthquakes  in  the 
California-Nevada  region.  Wyss  and  Brune  (1968)  found  that  earthquakes  in 
Nevada  often  generated  smaller  amplitude  surface  waves  for  a  given  Mi.  than  do 
earthquakes  along  the  San  Andreas  fault  or  along  the  northern  California  coast. 
In  a  later  study.  Wyss  and  Brune  (1971)  noted  large  differences  in  the  ratio  of 
Gutenburg  energy  Eq  determined  from  Mi.  with  seismic  moment  Mo  determined 
from  surface  waves,  for  earthquakes  in  the  California  region.  In  particular,  they 
found  high  values  of  apparent  stress  {/jlEq  / M0  )  for  earthquakes  in  the  region 
north  of  Bishop,  California.  At  that  time  it  was  not  possible  to  assess  the  effects 
of  variations  in  source  depth  and  source  mechanism. 

This  study  is  concerned  with  the  relative  excitation  of  seismic  energy  in  the 
band  0.05-10.0  Hz  from  earthquakes  in  the  1980  Mammoth  Lakes,  California, 
sequence.  These  earthquakes  have  been  studied  more  thoroughly  than  perhaps 
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any  other  earthquake  sequence  in  the  western  United  States.  A  large  number  of 
broad-band  digital  and  strong-motion  recordings  have  been  collected  in  the  epi- 
central  region,  and  a  dense  network  of  high-gain  seismographs  provide  data  for 
accurate  hypocentral  coordinates  and  first  motions  for  most  events.  The  epi- 
central  region  is  surrounded  by  a  network  of  long-period  seismographs,  and  is  at 
a  distance  of  20° -30°  from  the  high  gain  stations  of  the  Canadian  Network. 
Finally,  in  the  context  of  seismic  discrimination,  the  events  are  in  a  similar  tec¬ 
tonic  environment  to  the  Nevada  Test  Site,  approximately  150  km  to  the 
southeast  Thus  the  effects  of  structure  on  seismic  wave  propagation,  especially 
to  the  east,  should  be  similar  for  both  earthquakes  at  Mammoth  Lakes  and 
explosions  at  the  Nevada  Test  Site. 


Data 

This  study  is  based  on  seismograms  of  the  Mammoth  Lakes  earthquakes 
shown  in  Figure  1.  These  events  comprise  only  a  small  fraction  of  the  total 
number  of  events  which  have  occurred  in  the  Mammoth  Lakes  area  since  the 
earthquake  sequence  began  in  October,  1978.  Those  selected  for  study  consist 
ones  for  which  both  strong-motion  or  digital  local  recordings  and  long-period 
regional  recordings  exist.  All  events  have  been  located  using  arrival  times  from 
the  local  stations  shown  in  Figure  1,  with  additional  readings  from  local  stations 
operated  by  the  University  of  Nevada  and  the  California  Institute  of  Technology. 
We  have  obtained  first-motion  focal  mechanisms  for  a  number  of  the  events  from 
Vetter  (personal  communications,  1984). 

Seismograms  used  in  this  study  are  from  the  seismograph  stations  shown  in 
Figure  1.  The  SMA-1  strong-motion  accelerographs  were  installed  by  the  Califor¬ 
nia  Division  of  Mines  and  Geology  (Turpin,  1980),  or  by  L  ■  University  of  Southern 
California  (Moslem  et  at. ,  1983).  Digital  stations  described  by  Archuleta  et  at. 
(1982),  consist  of  either  force-balance  accelerometers  or  velocity  transducers 
and  Sprengnether  DR-100  recorders.  These  stations  are  denoted  by  triangles  in 
Figure  1.  For  events  along  the  southern  boundary  of  the  caldera,  the  SMA-1  and 
digital  seismographs  cover  a  larger  fraction  of  the  focal  sphere  than  many  previ¬ 
ous  studies  of  a  similar  nature.  The  radiation  pattern  of  events  south  of  the  cal¬ 
dera  is  not  as  well  sampled,  however. 

Figure  1  also  shows  the  locations  of  the  regional  long-period  stations 
(denoted  by  squares)  from  which  seismograms  were  used  to  determine  surface 
wave  magnitudes.  These  stations  completely  surround  the  Mammoth  Lakes 
region  and  are  all  WWSSN  stations  with  the  exception  of  the  California  Institute 
of  Technology  station  at  Pasadena  (PAS).  Readings  from  the  Canadian  network 
stations  (denoted  by  stars)  were  used  to  determine  body  wave  magnitude. 


Magnitude  Analysis 

Local  magnitude  ( ML  )  for  the  main  shocks  and  larger  aftershocks  are 
reported  from  Wood-Anerson  amplitude  measurements  at  the  University  of 
California-Berkeley  and  the  California  Institute  of  Technology.  Mi  for  the 
smaller  events  are  taken  from  Chavez  and  Priestley  (1985). 

The  Earthquake  Data  Reports  of  the  U.S.  Geological  Survey  report  peak-to- 
peak  amplitude,  period,  and  individual  station  body-wave  magnitude  (mb). 
From  this  listing  we  have  tabulated  the  body-wave  magnitude  reported  by  Cana¬ 
dian  Network  stations  at  distances  greater  than  20°  from  Mammoth  Lakes,  and 
from  these  values  we  have  determined  an  average  mb  value  for  the  37  Mammoth 
Lakes  earthquakes  studied. 
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We  have  determined  the  surface  wave  magnitude  using  the  definition  of 
Marshall  and  Basham  (1972).  This  value  of  Ma  is  defined  by  the  equation 

Ms=logA  +B  (A)  +  P(T)  (1) 

where  A  is  one-half  the  peak-to-peak  ground  amplitude  of  the  maximum  Rayleigh 
wave  with  a  period  T  at  a  distance  A.  5  (A)  is  a  correction  for  the  average 
effects  of  attenuation,  scattering,  geometrical  spreading  and  refraction;  P(T)  is 
a  correction  factor  for  dispersion  which  allows  measuring  the  surface  wave 
amplitude  over  a  range  of  periods.  Marshall  and  Basham  (1972)  have  tabulated 
P(T)  for  continental  North  America  for  periods  ranging  from  10  to  40  seconds. 
In  some  cases  the  maximum  surface  wave  amplitude  observed  for  the  Mammoth 
Lakes  earthquakes  was  in  the  range  8  to  10  seconds,  however  since  P(T)  is  a 
smoothly  varying  function,  we  have  extrapolated  P(T)  to  8  seconds  period. 
Marshall  and  Basham  have  also  provided  for  correcting  A/s  for  source  depth  by 
comparing  Ms  determined  for  an  event  at  different  periods.  However,  most  of 
the  surface  wave  observations  made  here  consist  of  a  single  Airy  phase,  and  thus 
such  a  comparison  could  not  be  made. 

Long-period  seismograms  from  all  western  United  States  WWSSN  stations 
and  the  long-period  station  at  Pasadena  were  examined.  Equ  (l)  was  used  to 
determine  the  Ma  from  eqn  (l)  for  each  station.  From  these  values,  we  deter¬ 
mined  an  average  Ma  for  each  event. 

Figure  2  compares  the  mb~Ms  to  the  North  American  earthquake  and 
explosion  data  of  Marshall  and  Basham  (1972).  Magnitudes  in  the  latter  study 
were  determined  from  seismograms  recorded  on  Canadian  Network  stations.  In 
this  study  Canadian  Network  data  were  used  to  determine  Tnb,  while  data  from 
long -period  stations  in  the  western  U.S.  were  used  to  determine  Ms.  This  could 
lead  to  a  systematic  difference  between  the  Ms  determined  by  Marshall  and  Bas¬ 
ham.  and  the  Afs  determined  for  the  Mammoth  Lakes  earthquakes.  This  does  not 
seem  to  be  a  strong  effect,  since  the  paths  overlap  to  a  large  extent,  and  the 
data  sets  in  Figure  2  overlap  as  well.  Neither  the  Marshall-Basham  data  nor  the 
Mammoth  Lakes  earthquakes  data  have  been  corrected  for  depth  variations.  All 
Mammoth  Lakes  events  fall  within  the  earthquake  population  and  discriminate 
from  explosions,  although  some  events  are  well-offset  from  the  earthquake  mean 
towards  the  explosion  population. 


Body- wave  estimates  of  source  parameters 

The  seismograms  analyzed  display  varying  degrees  of  complexity,  ranging 
from  simple  pulses  to  wave  trains  several  seconds  long.  This  variation  is  possibly 
a  consequence  of  the  position  of  the  station  on  the  focal  sphere.  For  instance, 
seismograms  showing  large,  simple  S-waves  often  have  relatively  small  P-wave 
amplitudes.  In  this  case,  if  the  S-wave  amplitude  is  large,  it  will  exceed  the  scat¬ 
tered  energy  and  the  resulting  seismogram  will  have  a  simpler  character.  In 
some  cases,  it  is  apparent  from  the  long  period  oscillations  following  the  S-wave, 
that  earth  structure  contributes  to  the  complexity  of  the  seismograms.  Dis¬ 
placement  amplitude  specta  were  computed  for  these  locally  recorded  body- 
wave  data.  The  SMA-1  analog  seismograms  were  digitized  at  a  sample  interval  of 
0  025  sec;  the  digital  recordings  supplied  from  the  U.S.  Geological  Survey  were 
digitized  at  an  interval  of  0.05  sec.  Displacement  amplitude  spectra  were  com¬ 
puted  by  Fourier  transforming  the  S-wave  arrival,  correcting  for  instrument 
response,  and  dividing  the  acceleration  spectra  by  cj-2  and  the  velocity  spectra 
by  cj-'  to  give  displacement  spectra.  Time  windows  were  chosen  to  be  long 
enough  to  include  the  entire  S-wave  arrival,  but  short  enough  to  minimize  the 
inclusion  of  surface  wave  energy.  The  spectra  have  been  normalized  to  a 
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distance  of  10  km.  Rather  than  rotate  the  data  into  radial  and  transverse  com¬ 
ponents,  we  took  the  vector  sum  of  the  two  horizontal  component  spectra. 

The  spectra  were  interpreted  in  terms  of  a  flat,  low-frequency  asymptote, 
do,  a  high  frequency  falloff  with  a  slope  in  the  range  cj-1  to  cj~3,  and  the  corner 
frequency  fe  at  which  the  two  trends  intersect.  In  almost  all  cases,  there  is  a 
second  corner  in  the  data  at  frequencies  between  10  and  16  Hz.  Below  this  fre¬ 
quency  the  spectra  drop  off  as  cj~*  to  cj-6.  This  frequency  may  represent  the 
/max  of  Hanks  (1982).  The  asymptotes  were  fit  by  eye  using  an  interactive  com¬ 
puter  routine.  This  method  of  fitting  is  subjective  and  a  more  rigorous  tech¬ 
nique  could  be  adopted  such  as  fitting  the  asymptotes  by  a  least-squares  type 
procedure.  However,  the  simple  technique  used  here  is  felt  to  be  adequate 
because  its  errors  were  smaller  than  other  known  errors  involved  in  estimating 
Q„  and  fc.  Figure  3  is  an  exampLe  of  this  estimation  of  the  fits  for  one  seismo¬ 
gram. 


The  far-field  spectral  parameters  were  then  related  to  the  source  parame¬ 
ter  seismic  moment  M0  using  the  relationship  given  by  Brune  (1970,  1971) 


_  4np/S3(l0r 

0  ~ 
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where  p  is  the  density  (p  =  2, 8  <7  cm-3),  &  is  the  shear-wave  velocity 

(/?  =  3.5jtms"'),  r  is  the  hypocentral  distance,  is  the  RMS  average  of  the 
radiation  pattern  (taken  ha”e  as  0.6),  and  fc  is  a  correction  factor  for  the 
amplification  of  SH  waves  u  the  free  surface  (taken  as  2).  The  correction  of  the 
amplitude  of  the  S-«u.  on  being  reflected  at  the  free  surface  assumes  that  the 
waves  are  horizon'..?'  polarized  (SH),  however  since  we  have  not  rotated  the 
data  into  transverse  and  radial  components,  the  data  are  in  general  some  mix¬ 
ture  of  SH  a  id  SV  motion.  Abe  (1974)  discusses  experimental  studies  of 
amplificati:  u  on  reflection  and  concludes  that  using  fc  =  2  applies  approximately 
to  SH-wa'  e  signals  containing  some  SV  motion.  Helmberger  and  Malone  (1975) 
discuss  problems  due  to  multiple  arrivals  caused  by  reflections  in  deducing 
source  parameters  from  S-waves  recorded  at  epicentral  distances  much  greater 
than  the  source  depth.  The  station-source  configuration  shown  in  Figure  1  indi¬ 
cates  that  this  should  not  be  a  problem  in  this  study. 

A  source  of  error  inherent  in  the  estimate  of  fl0,  and  consequently  in  the 
estimates  of  M0,  comes  from  the  assumption  of  a  single,  average  correction  for 
radiation  pattern.  If  sufficient  instruments  were  available,  or  if  the  source 
parameters  including  direction  and  velocity  of  rupture  were  known,  this  error 
could  be  reduced.  A  value  of  0.6  was  used  for  R^  (Thatcher  and  Hanks,  1973). 
At  some  azimuths,  the  amplitude  could  be  1.67  times  larger  than  the  amplitude 
determined  assuming  this  correction,  while  near  a  node  in  the  radiation  pattern 
the  amplitude  could  be  quite  small.  Furthermore,  if  the  source  dimension  is 
comparable  to  the  wavelength  involved,  source  propagation  might  distort  the 
radiation  pattern  by  focusing  or  defocusing  energy,  thereby  changing  the 
apparent  corner  frequency  for  different  directions  and  possibly  reducing  the 
calculated  moment. 

No  Q  correction  was  made  in  determining  the  spectra.  Archuleta  et  at. 
(1982)  found  that  whole-path  Q  had  a  small  effect  on  the  measurement  of  Qc  for 
Mammoth  Lakes  events. 

Errors  in  f!0  (  and  consequently  in  M0  )  were  estimated  by  calculating  an 
average  spectral  value  and  a  standard  deviation  where  more  than  one  recording 
was  available  for  an  event.  The  average  values  of  fi0  were  determined  by  the  for¬ 
mula 
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where  Ho.  is  the  long-period  spectral  level  at  the  ith  station  ,  Rx  is  the  distance 
to  the  ith  station,  and  NS  is  the  number  of  stations  recording  the  event.  The 
standard  deviation  of  the  log  average  value  of  fi0  was  determined  by  the  formula 
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A  multiplicative  error  factor  is  then  determined  from 

EQ o  =  antilog  js,  d,  ( log  <  G0  > ) 


(5) 


As  pointed  out  by  Archuleta  et  al.  (1982),  when  calculating  average  values  of  Cl0 
(or  Mq),  it  is  necessary  to  compute  the  averages  using  the  equations  in  the 
above  forms  to  give  equal  weight  to  each  observation.  If  a  simple  arithmetical 
average  value  is  determined  it  will  be  biased  toward  larger  values,  since  the 
errors  associated  with  n„  are  lognormally  distributed. 


Determination  of  Surface  Wave  Moment 

Tucker  and  Brune  (1977)  used  surface  wave  spectra  of  the  San  Fernando 
aftershocks  to  obtain  estimates  of  seismic  moments  which  were  independent  of 
those  indicated  by  the  S-wave  spectra,  and  which  sampled  a  lower  frequency 
portion  of  the  spectrum.  For  a  number  of  large  eventsmoments  determined  at 
the  surface  wave  periods  indicated  that  the  low-frequency  part  of  the  spectrum 
is  roughly  proportional  to  cj-1.  Reexamination  of  the  near-source  spectra  for 
these  events  were  in  fact,  found  to  be  rising  towards  lower  frequency  in  contrast 
to  the  spectra  for  almost  all  other  events.  They  made  a  careful  estimate  of  the 
errors  involved  in  estimating  the  seismic  moment  from  the  surface  wave  spec¬ 
tra.  The  errors  included  procedural  errors,  path  structure  and  attenuation, 
lateral  refraction  (multipathing),  source  finiteness,  and  uncertainties  in  fault 
parameters  (depth,  dip,  strike,  and  slip  vector)  in  correcting  for  the  radiation 
pattern.  They  concluded  that  these  effects  could  lead  to  a  maximum  error  in 
the  final  result  of  no  more  than  a  factor  of  4. 

In  this  study  we  also  determined  the  seismic  moment  from  surface  waves 
for  comparison  with  seismic  moments  determined  from  the  S-wave  spectra,  how¬ 
ever  we  used  a  different  procedure  which  allowed  us  to  determine  the  seismic 
moment  from  a  large  number  of  events  without  the  time-consuming  digitization 
of  the  surface  wave  seismograms.  Surface  wave  synthetics  were  generated  by  a 
computer  routine  (written  by  Guy  Masters)  for  an  earthquake  with  M0  =  1021 
dyne-cm  at  a  number  of  regional  seismographs  from  which  data  were  used  to 
determine  the  surface-wave  magnitudes.  From  the  synthetic  surface  wave 
seismograms,  we  determined  the  surface  wave  magnitude  in  the  same  manner 
as  discussed  above  for  the  Mammoth  Lakes  earthquakes,  and  thus  determined  a 
linear  relationship  between  the  surface  wave  magnitude  and  the  logarithum  of 
the  seismic  moment. 

Factors  affecting  the  calculation  of  the  synthetic  surface  wave  seismograms 
are  the  earth  velocity  and  attenuation  structure  and  the  fault  parameters  of  the 
individual  earthquakes.  Since  the  events  are  located  at  the  western  edge  of  the 
Great  Basin,  while  most  of  the  seismograph  stations  are  located  in  western  North 
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America,  we  have  used  the  velocity  structure  determined  for  the  Great  Basin 
from  surface  wave  dispersion  measurements  (Priestley  and  Brune,  1979).  The 
attenuation  structure  for  the  Great  3asin  by  Patton  and  Taylor  (1984)  was  used, 
with  a  slight  lowering  of  Qp  from  that  of  Patton  and  Taylor  in  the  surficial  layer. 
Since  we  do  not  have  fault-plane  solutions  nor  do  we  know  the  focal  depth  for  the 
majority  of  events  we  have  studied,  we  determined  an  average  relationship  for 
the  Mammoth  Lakes  area  in  the  following  manner.  Vetter  and  Ryall  (1983)  have 
determined  fault  plane  solutions  for  a  large  number  of  events  in  the  Mammoth 
Lakes  area,  and  found  that  the  mechanisms  can  generally  be  divided  into  two 
groups.  The  first  group  consists  of  earthquakes  with  strike-slip  mechanisms 
(average  plane  of  strike  N14°E,  dip  80°ESE,  slip  N16°E  or  strike  \73°W,  dip 
80°SS'W,  slip  N76°E)  which  are  most  commonly  at  depths  less  than  9  km.  The 
second  group  consists  of  earthquakes  with  oblique  or  normal  mechanisms  (aver¬ 
age  planes  of  strike  N03°E,  dip  64°E,  slip  N40°E  or  strike  NT50°W,  dip  380SW,  and 
slip  \86°W)  usually  at  depths  greater  than  9  km.  Synthetic  surface  wave  seismo¬ 
grams  were  calculated  for  both  mechanisms  for  2,  5,  8,  11  and  14  km  depth. 
Surface-wave  magnitudes  were  determined  from  the  10  records  at  each  station 
and  an  average  conversion  factor  determined  for  both  focal  mechanisms  and  all 
focal  depths.  The  resulting  relationship  between  surface  wave  magnitude  Ms 
and  seismic  moment  M 0  for  the  Mammoth  Lakes  area  is  given  by  the  equation 

log  Ma  =  Ms  +  19.40(±0.06)  .  (5) 

To  verify  the  validity  of  our  moment-surface  wave  magnitude  relationship, 
we  can  compare  equation  (5)  with  other  similar  relations  established  for  the 
western  United  States.  Wyss  and  Brune  (1968)  originally  established  a  moment- 
local  magnitude  relationship  for  events  less  than  magnitude  6  in  the  western 
U.5.  Seismic  moments  were  determined  for  272  events  which  averaged  over  a 
wide  range  of  tectonic  regions  (California,  Nevada,  Arizona,  Utah  and  Baja  Cali¬ 
fornia).  These  were  found  to  be  besf  fit  by  the  curve  given  by  the  equation 

log  M0  =  1.7  Mi  +  15. 1  3  <  Ml  <  6  (7) 

Using  the  Ms  -  Mi  relationship  g:ven  by  Wyss  and  Brune  (1968),  this  becomes  in 
terms  of  MS(C  )  (the  surface  wave  magnitude  defined  by  Gutenberg  (1945)) 

log  Mq  =  MS(C)  +  19.2  3  <  Ms  <  6  .  (9) 

This  definition  can  be  compared  with  the  relationship  derived  directly  from 
Gutenberg's  (1945)  definition  of  surface  wave  magnitude.  According  to  that 
definition,  a  magnitude  6  earthquake  produces  a  far-field  displacement  of 
lOOjxm  at  distance  22°  for  surface  waves  of  period  20  seconds.  In  terms  of 
moment  Gutenburg's  relationship  becomes 

log  M0  =  MS(G)  +  19.3  (9) 

This  point  thus  represents  average  of  the  numerous  observations  on  which  the 
surface  wave  magnitude  was  based.  As  pointed  out  by  Richter  (1958),  the  scale 
was  adjusted  to  agree  with  the  local  magnitude  values  of  6  to  7.  Tucker  and 
Brune  (1977)  did  not  explicitly  give  a  relationship  between  surface  wave  magni¬ 
tude  and  seismic  moment,  however  from  the  values  given  in  their  Tables  1  and  3, 
the  average  relationship  is 


log  M 0  =  Ms  +  19.3 


(10) 
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This  relationship,  determined  using  a  Gutenburg  continental  earth  model  and 
the  Marshall-Basham  definition  of  surface  wave  magnitude,  is  identical  to  the 
relationship  based  on  the  Gutenburg  definition  of  surface  wave  magnitude.  Com¬ 
paring  these  various  relationships  between  Afs  and  M0  we  conclude  that  the  use 
of  the  Great  Basin  velocity  (Priestley  and  Brune,  1978)  and  attenuation  (Patton 
and  Taylor,  1984)  model  has  increased  the  estimated  moment  of  the  events  by 
26%  over  that  determined  for  the  Gutenburg  continental  model. 

Equation  (6)  aLlows  us  to  rapidLy  determine  the  seismic  moment  Ma  from 
the  surface  wave  magnitude  observation.  The  resulting  values  of  seismic 
moment  were  converted  to  an  equivalent  value  of  Dc  by  solving  equation  (2)  for 
Go  and  substituting  the  seismic  moment  derived  from  the  surface  wave  observa¬ 
tions.  This  gives  values  of  G0  at  frequencies  of  the  surface  waves  for  comparison 
with  the  spectra  obtained  from  the  body-wave  record  on  the  accelerometers. 
The  values  of  seismic  moment  are  compared  in  Tapbe  4. 


Comparison  of  short-period  and  long-period  estimates  of  seismic  moment 

Figure  4  is  a  plot  of  the  seismic  moment  determined  from  the  spectral 
analysis  of  the  strong-motion  and  digital  data  recorded  at  near  distances  (A  <  20 
km),  and  the  seismic  moment  determined  from  the  surface  waves  recorded  at 
larger  distances  (A  >  375  km).  The  line  is  of  unit  slope  and  zero  intercept,  and  is 
for  reference.  As  shown  in  this  Figure,  most  of  the  events  we  have  studied  show 
no  significant  differences  in  the  long-period  spectral  level  as  measured  near  fre¬ 
quencies  near  1  Hz  and  as  measured  near  frequencies  of  10  Hz.  1.0  Hz.  In  com¬ 
paring  short  period  versus  long  period  excitation  of  events  in  the  California 
region,  Wyss  and  Brune  (1971)  found  that  events  occurring  to  the  north  of 
Bishop,  California  indicated  higher  apparent  stress  than  events  along  the  San 
Andreas  fault.  In  a  number  of  recent  studies  in  southern  California  (Tucker  and 
Brune.  1977;  Hartzell  and  Brune,  1977,  1979),  discrepancies  have  been  note 
between  the  between  long  and  short  period  determinations  of  seismic  moments. 
Such  discrepancies  may  arise  due  to  partial  stress  drop  (Brune,  1970)  or  to  aft¬ 
erslip  on  the  fault  (Hartzell  and  Brune,  1979).  In  the  latter  case,  the  body  wave 
moment,  and  presumably  mb  and  ML,  are  the  result  of  the  breaking  of  an  asper¬ 
ity,  while  the  overall  moment,  and  presumably  Ms,  is  the  result  of  the  overall 
faulting  process.  We  have  found  such  events  to  be  uncommon  in  the  Mammoth 
Lakes  sequence,  as  is  shown  in  Figure  4.  This  may  arise  due  to  a  higher  overall 
state  of  stress  in  the  crust  in  the  vicinity  of  Mammoth  Lakes. 


Comparison  of  magnitude  estimates  with  source  spectra 

Six  events  were  chosen  for  special  study.  These  are  events  which  cover  the 
range  of  depths  and  focal  mechanisms  observed  at  Mammoth  Lakes,  and  for 
which  we  have  local  data  for  determination  of  the  spectra,  plus  determinations 
of  ML  ,  m6  ,  and  Ms  . 

In  order  to  compare  the  similarities  and  differences  in  the  magnitude  esti¬ 
mates  for  these  events,  it  is  important  to  consider  the  frequency  associated  with 
each  measurement.  The  magnitude  values  were  related  to  the  spectra  in  the  fol¬ 
lowing  way.  The  301481901  event  was  chosen  as  a  reference  since:  (l)  both  the 
spectrum  and  magnitude  estimates  were  well-determined,  (2)  this  was  one  of  the 
smaller  events  and  thus  any  effects  of  fault  finiteness  are  minimized;  (3)  the 
event  is  shallow;  and  (4)  the  frequencies  at  which  Mi,  mh,  and  Ms  were  deter¬ 
mined  at  are  all  lower  than  the  corner  frequency,  i.e.,  all  measurements  were 
from  the  flat  portion  of  the  spectrum.  The  average  spectrum  for  this  event  is 
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shown  in  Figure  5a.  The  error  bars  at  the  corner  of  the  spectrum  denote  the 
range  in  0a  and  fc ,  and  the  high  frequency  slope  is  the  average  of  all  measure¬ 
ments.  Additive  constants  for  each  of  the  magnitudes  were  found  such  that  the 
magnitude  values  fall  on  the  spectrum  corresponding  to  their  appropriate  fre¬ 
quency.  The  frequency  for  m6  and  Ms  is  the  mean  frequency  of  the  magnitude 
measurements  and  error  bars  denote  the  spread  in  those  measurements.  Mi 
was  assumed  to  be  measured  at  a  period  of  0. B=0.2  sec.,  the  free  period  of  a 
Wood-Anderson  seismometer 

The  magnitude  measurements  were  made  at  similar  frequencies  for  all  six 
earthquakes.  The  most  significant  deviation  is  for  m b  where  the  period  for  the 
shallower  events  is  0.3  =  0. 1  sec  and  the  period  for  the  deepest  event,  301471224, 
is  1.3  =  0. 1  sec.  Because  all  paths  to  the  regional  and  teieseismic  stations  were 
essentially  identical,  and  because  of  the  similarity  in  the  frequency  of  each  of 
the  magnitude  measurements,  we  felt  justified  in  applying  the  additive  constants 
determined  for  the  301481901  event  to  the  other  five  events. 

The  301490516  event  is  similar  in  magnitude  and  depth  to  the  reference 
event.  Both  events  have  a  similar  corner  frequency  and  it  is  not  surprising  that 
Mi  and  mb  fall  below  the  corner  frequency  and  on  the  flat  portion  of  the  spec¬ 
trum.  The  Ms  value  falls  above  the  long-period  level,  but  this  is  not  significant 
considering  the  error  bars  on  Ms  and  fl„. 

The  event  301471357  is  almost  one  magnitude  unit  larger  than  the  refer¬ 
ence  event  801431901,  but  is  of  similar  depth.  The  Ms  value  for  this  event  falls 
close  to  the  long-period  level  while  both  mb  and  Mi  were  measured  at  frequen¬ 
cies  within  the  error  bars  fc.  When  these  values  are  normalized  to  the  refer¬ 
ence  event,  both  fall  below  the  long  period  level  and  near  the  high  frequency 
slope. 

Comparison  of  these  three  events,  801481901  and  301490516  which  are  of 
similar  magnitude  and  depth,  and  801471857  which  is  of  similar  depth  but 
approximately  one  magnitude  unit  larger  justifies  our  rationale  of  relating  the 
magnitude  measurements  to  the  spectrum. 

The  event  occurring  at  801462035  has  the  largest  surface  wave  magnitude 
of  the  six  event  discussed  in  this  section.  It  is  approximately  3  km  deep,  near 
the  transition  zone  in  focal  mechanisms  found  by  Vetter  (1984).  We  have  only 
one  strong-motion  record  for  determining  the  source  spectrum,  hence  there  are 
no  error  bars  on  the  corner  frequency  or  long-period  level.  mb  is  determined  at 
a  frequency  lower  than  the  corner  frequency,  and  falls  approximately  a  factor  of 
two  below  the  fi0  level.  Mi  is  above  the  corner  frequency  and  not  far  removed 
from  the  high  frequency  asympote.  The  Ms  value,  however,  is  a  factor  of  10 
greater  than  the  long-period  asymptote.  This  event  is  approximately  6  km 
deeper  than  those  discussed  above  which  would  somewhat  decrease  the  short 
period  Rayleigh  wave  excitation.  The  mean  period  of  the  Rayleigh  wave  Airy 
phase  is  slightly  longer  for  this  event  (I2.0±2.5sec)  compared  to  the  shallower 
events  ( 10.0±2.0sec).  Such  discrepancies  between  long  and  short  period 
moments  have  been  observed  in  a  number  of  similar  studies  in  southern  Califor¬ 
nia  (Tucker  and  Brune,  1977;  Hartzell  and  Brune,  1977,  1979).  Such  a 
discrepancy  may  arise  due  to  partial  stress  drop  (Brune,  1970)  or  to  afterslip  on 
the  fault  (Hartzell  and  Brune,  1979).  In  the  latter  case,  the  body  wave  moment, 
and  presumably  mb  and  Mi ,  are  the  result  of  the  breaking  of  an  asperity,  while 
tbe  overall  moment,  and  presumably  Ms ,  is  the  result  of  the  overall  faulting  pro¬ 
cess.  We  have  found  such  events  to  be  uncommon  in  the  Mammoth  Lakes 
sequence,  as  is  shown  in  Figure  4. 

The  event  301471224  occurring  at  12.5  km,  is  the  deepest  event  studied,  as 
well  as  one  of  the  deepest  events  of  the  Mammoth  Lakes  sequence.  For  this 
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event,  Mi  falls  on  the  high  frequency  asymptote  above  the  corner  frequency, 
while  7nb  is  at  a  frequency  below  the  corner  and  approximately  a  factor  of  five 
beneath  the  long  period  level.  The  Ms  error  bars  intersect  the  long  period  level 
although  the  average  value  is  slightly  below,  possibly  due  to  the  increased 
source  depth. 

The  final  event  examined  is  that  occurring  at  801520050.  This  is  the  smal¬ 
lest  event  of  the  six  studied  and  one  of  those  in  Figure  2  which  is  offset  from  the 
earthquake  population  towards  the  explosion  population.  The  depth  of  this 
event  is  similar  to  that  of  the  reference  event,  so  the  low  value  of  Ms  cannot  be 
attributed  to  source  depth.  Both  mb  and  Mi  fall  within  the  corner  frequency 
error  bars.  Two  features  of  the  spectrum  of  this  event  distinguish  it  from  the 
other  five  events  examined.  First,  the  high-frequency  slope  is  nearer  to  :j_1 
compared  to  a  high  frequency  slope  of  approximately  u-2  for  the  other  events. 
Second,  in  all  of  the  individual  spectra  for  801520058,  there  is  a  sharp  break  in 
slope  at  about  10  Hz.  The  other  earthquakes  also  show  a  second  break  in  slope 
between  10  and  15  Hz  which  may  represent  the  same  feature.  This  may  be  more 
pronounced  for  the  801520058  event  due  to  the  shallower  slope  between  1  and  10 
Hz.  Other  than  these  features,  there  is  nothing  to  distinguish  this  event  from 
those  which  plot  along  the  earthquake  trend  in  Figure  2.  All  of  the  six  events 
with  the  exception  of  801520058  have  stress  drops  of  the  order  100  bars. 
801520058  has  a  stress  drop  of  15  bars. 

The  moment  and  corner  frequencies  for  the  801471857,  801481901  and 
801490516  events  are  all  in  approximate  agreement  with  the  cj-2  model  for  the 
spectra  of  earthquakes  (Aki,  1967;  Brune  1970,  1971).  The  corner  frequency  of 
the  801471857  earthquake  is  about  a  factor  of  2  lower  than  801481901  and 
801490516,  whereas  the  corresponding  long-period  amplitudes  of  801471857  are 
approximately  10  times  higher.  Thus  the  corner  frequencies  and  moment  of  the 
801471857,  801181901,  and  801490516  are  in  agreement  with  the  model  for 
events  of  approximately  the  same  stress  drop,  but  with  source  dimensions 
different  by  a  factor  of  2  (Hanks,  1979). 


Conclusions 

In  comparing  short  period  versus  long  period  excitation  of  events  in  the  Cal¬ 
ifornia  region,  Wyss  and  Brune  (1971)  found  that  events  occurring  to  the  north  of 
Bishop,  California  indicated  higher  apparent  stress  than  events  along  the  San 
Andreas  fault.  In  this  study  we  have  looked  in  detail  at  the  relative  excitation  of 
short  and  long  period  energy  in  this  region  from  aftershocks  of  the  Mammoth 
Lakes  sequence. 

We  have  compared  the  body  wave  magnitude.  rnb ,  and  surface  wave  magni¬ 
tude,  Ms ,  and  find  that  Mammoth  Lakes  events  with  few  exceptions  discriminate 
from  nuclear  explosions  occurring  at  the  Nevada  Test  Site  approximately  150 
km  to  the  southeast. 

We  have  used  local  strong-motion  and  broad-band  digital  recordings  to  com¬ 
pute  the  source  spectra  of  37  earthquakes.  For  events  in  common  with 
Archuleta  et  a L  (1982)  our  moments  are  in  agreement.  We  have  developed  a 
moment-surface  wave  magnitude  relationship  which  allows  us  to  rapidly  deter¬ 
mine  the  seismic  moment  at  surface  wave  periods  from  time  domain  amplitude 
measurements  on  the  long-period  seismogram.  These  determinations  of 
moment  in  almost  all  cases  agree  closely  with  the  moments  determined  from 
the  spectra.  Similar  studies  in  southern  California  show  long-period  moments 
exceeding  short-period  moments  by  a  factor  of  4  to  10. 
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Figure  Captions 

Figure  1.  Location  maps  showing  the  events  studied  and  the  stations  from  which 
seismograms  were  analysed.  In  the  lower  right-hand  map,  the  solid  circles 
denote  locations  of  events,  triangles  denote  location  of  strong-motion 


-  12  - 


accelerographs,  and  digital  seismographs.  The  upper  left-hand  map  shows 
the  locations  of  the  long-period  seismographs  used  for  the  surface  wave 
analysis  (squares),  and  Canadian  network  short-period  stations  used  to 
determine  body-wave  magnitude. 


Figure  2.  Comparison  of  the  body-wave  magnitudes  surface-wave  magnitude 
The  solid  dot  (earthquakes)  and  pluses  (explosions)  are  taken  from 
Marshall  and  Basham  (1972);  the  solid  squares  are  the  Mammoth  Lakes 
events  from  this  study.  All  body-wave  magnitudes  are  from  measurements 
from  Canadian  network  seismographs.  Surface  wave  magnitudes  taken 
from  Marshall  and  Basham  are  from  Canadian  network  long -period  seismo¬ 
graphs,  surface  wave  magnitudes  for  the  Mammoth  Lakes  earthquakes  are 
from  long-period  WWSSN  stations  in  the  western  United  States. 


Figure  3.  Representative  analysis  of  short-period  digital  seismogram.  The 
bracket  over  the  seismogram  indicates  the  S-wave  window;  the  heavy  solid 
lines  indicate  the  long-period  and  high-frequency  asymptotes  fit  to  the 
spectrum. 


Figure  4.  Comparison  of  the  body-wave  and  surface-wave  determined  moments 
for  the  events  studied.  The  line  is  for  reference  and  is  of  unit  slope  and 
zero  intercept. 


Figure  5.  Comparison  of  the  source  spectra  and  the  determined  values  of  the 
local  magnitude,  body-wave  magnitude,  and  the  surface  wave  magnitude  at 
threir  appropriate  frequencies.  The  error  bars  assigned  to  each  spectral 
value  represent  the  estimated  standard  error  of  the  measurement. 
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THE  Ml  SCALE  IN  THE  GREAT  BASIN  AND  M0  vs  ML 
RELATIONSHIPS  FOR  THE  1980  MAMMOTH  LAKES, 
CALIFORNIA  EARTHQUAKE  SEQUENCE 


Davn.dE.  Chavez  and  Keith  F  Priestley 


ABSTRACT 


We  have  converted  471  seismograms  from  103  earthquakes  in  the 
vicimty  of  Mammoth  Lakes,  California  into  equivalent  Wood- 
Anderson  seismograms  and  estimated  local  magnitude  Ml  following 
the  original  definition  made  by  Richter,  We  found  that  the  distance 
correction  terms  given  by  Richter’s  attenuation  curve  yield  sys¬ 
tematically  larger  Ml  estimates  for  distances  less  than  20  km.  Ml 
determined  at  near  distances  (  AclOkm  )  using  Richter's  logA0 
values  are  as  much  as  1  magnitude  unit  greater  than  M;t  deter¬ 
mined  at  regional  distances  {  A >50  km  )  for  the  same  event.  This 
result  is  similar  to  those  of  recent  studies  utilizing  southern  Cali¬ 
fornia  and  northern  Baja  California  eartnquakes.  Seismograms 
-acorded  at  greater  distances  m  the  Great  Basin  (100  km  <A<  300 
.<m)  show  more  attenuation  (about  0.3  magnitude  units)  than  that 
predicted  by  Richter’s  curve.  We  have  determined  corrections  to 
the  attenuation  curve  which  removes  the  distance  dependence 
observed  for  our  data.  We  also  find  a  near  source  local  magnitude 
saturation  for  events  with  Ml  ^6.0  similar  to  earlier  studies  with 
California  earthquakes  and  explained  as  finite  source  size  effects. 
Using  the  revised  magnitude  scale  and  seismic  moments  estimated 
from  spectral  analysis  we  found  that  our  data  were  well  fit  by  the 
straight  line 

log  M0  =  (  1.22-0. C4)A/j,  -i-  ( 16.93±0. 15) 

for  1.5  §  d/j,  <  6.3.  Comparison  of  our  data  to  predicted  A/0  versus 
Ml  values  suggest  a  constant  stress  drop  of  approximatley  100 
bars  for  events  with  Ml  <5.0  and  increasingly  greater  stress  drop 
for  the  larger  events. 


Introduction 


Local  magnitude,  ML .  was  defined  by  Richter  (1925)  to  be 

Ml  =  log  A  -  logA„  (A) 

wnere  A  is  the  maximum  zero  to  peak  trace  amplitude  in  millimeters  recorded 
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on  a  standard  Wood-Anderson  torsion  seismograph  (static  magnification  =  2B00, 
natural  period  =  0.8  sec,  and  damping  factor  =  0.8)  at  an  epicentral  distance  A. 
The  logA„  (A)  function  removes  the  distance  dependence  of  the  amplitude,  incor¬ 
porating  geometrical  spreading,  change  in  wave  type,  scattering,  and  anelastic 
attenuation.  The  nature  of  these  last  three  will  depend  on  the  structure 
between  the  source  and  receiver  and  in  general  will  vary  between  regions. 
Richter  emperically  determined  the  shape  of  the  logA0  (A)  curve  for  southern 
California  and  calibrated  it  so  that  an  event  of  Mi  =0  would  have  a  peak  Wood- 
Anderson  trace  amplitude  of  one  micron  at  an  epicentral  distance  of  100  km. 
Initially  the  curve  was  defined  for  the  range  25 km  <  A  <  600 km,  however  later 
Gutenturg  and  Richter  (1942)  extended  it  to  zero  distance  using  data  from  low 
gain  torsion  seismographs. 

The  logA0  (A)  values  published  in  Richter  (1958,  p.  342)  have  been  used  in 
estimating  Mi  outside  of  the  southern  California  region  for  which  they  were 
developed.  As  pointed  out  in  a  number  of  studies,  this  implies  that  the  attenua¬ 
tion  function  is  the  same  as  in  southern  California.  This  will  not  always  be  the 
case,  however,  and  consequently  equal  numerical  Mi  values  for  earthquakes  in 
different  regions  will  not  neccesarily  imply  the  same  seismic  energy  release  at 
the  source.  This  is  of  interest  in  engineering  seismology  because  Mi  (and  hence 
log.40  (A)  )  plays  an  important  role  in  the  analysis  and  prediction  of  strong 
ground  motion. 

Few  Wood-Anderson  attenuation  functions  have  been  determined  for  areas 
outside  of  southern  California  primarily  because  of  the  lack  of  a  sufficient 
number  of  standard  seismographs  with  which  to  constrain  the  curve.  The 
increasing  availability  of  digital  seismograms  from  calibrated  stations  now  offers 
'■.he  possibility  of  doing  so  by  creating  "synthetic”  Wood-Anderson  seismograms. 
For  example.  Bakun  (1984a)  used  such  data  to  develop  a  logA0  (A)  curve  for  cen¬ 
tral  Cai  Forma.  He  found  that  the  original  southern  California  curve  was 
appropriate  except  for  some  modification  at  near  distances.  Also,  three  recent 
studies  in  the  southern  California  region  (Luco,  1982;  Jennings  and  Kamamori, 
1383;  Viunguia  and  Brune,  1984)  have  used  strong  motion  and/or  digital  record¬ 
ings  to  determine  Mi  at  near  instances.  Each  of  these  have  found  a  systematic 
trend  in  ML  with  distance,  indicating  that  the  !ogA0  attenuation  function  used  in 
calculating  Mi  is  not  correct  close  to  the  epicenter.  Moreover,  they  found  that 
the  nature  of  these  variations  with  distance  are  a  function  of  recorder  site 
characteristics  and  the  size  of  the  event  as  well.  For  instance,  all  reported  that 
for  large  events  {ML  >  8.0)  local  magnitude  was  underestimated  at  distances 
below  about  20  km  relative  to  that  determined  at  more  distant  stations. 
Munguia  and  Brune  also  found  that  for  smaller  events  (ML  §5.5)  local  magnitude 
was  increasingly  overestimated  as  epicentral  distance  grew  less.  In  all  cases,  no 
modificauons  to  Richter's  (1358)  curve  were  found  to  be  necessary  beyond  100 
km. 

In  this  paper  we  use  events  in  the  vicinity  of  Mammoth  Lakes,  California 
recorded  both  near  the  source  and  at  regional  distances  to  develop  a  log.4„  (A) 
curve  for  the  Great  Basin.  We  are  fortunate  in  that  ray  paths  from  Mammoth 
Lakes  eastward  travel  exclusively  in  the  Great  Basin  while  paths  westward  travel 
through  the  areas  where  Richter's  (1958)  curve  was  defined  or  later  shown  to  be 
applicable.  This  allows  us  to  directly  calibrate  our  attenuation  curve  to  that  for 
wnich  Mi  was  originally  defined. 

For  most  applications  outside  of  engineering,  the  seismic  moment  M0  is 
often  a  preferable  measure  of  earthquake  size  since  moment  is  directly  related 
to  specific  source  properties,  namely  the  fault  area  and  amount  of  slip.  More¬ 
over,  Mi  will  saturate  for  events  sufficiently  large  that  the  corner  frequency  is 


-  3  - 


below  the  Wood-Anderson  passband  while  M0  will  not.  Theoretical  source  models 
and  the  widespread  use  of  digital  recordings  have  made  possible  routine  estima¬ 
tion  of  M0  for  a  large  range  of  sizes.  Here,  we  develop  a  M0  versus  Mi  relation¬ 
ship  for  events  in  the  1980  Mammoth  Lakes  sequence  and  compare  it  to  similar 
relations  for  other  areas. 


Data  and  Analysis 

The  data  analyzed  in  this  study  consist  of  seismograms  of  -»  "thquakes  in 
the  vicinity  of  Mammoth  Lakes,  California,  primarily  events  the  i990  sequence 
(Urhammer  and  Ferguson,  1980).  The  data  at  distances  less  than  30  km  include 
velocity  and  acceleration  records  obtained  by  the  U.  S.  Geological  Survey 
(Archuleta,  et  al.,  1982),  and  strong  motion  accelerograms  from  instruments 
installed  by  the  California  Divisions  of  Mines  and  Geology  (Turpin,  1980)  and  by 
the  University  of  Southern  Caiiforma-California  Institute  of  Technology  (Moslem, 
et  al.,  1983).  For  distances  from  100  to  600  km  we  have  data  from  broad-band 
stations  operated  by  the  Lawrence  Livermore  National  Laboratory.  Similar 
broad-band  instruments  installed  by  the  University  of  Nevada  recorded  aft¬ 
ershocks  of  the  November  23,  1984  Mi  =5,7  earthquake  near  Bisnop,  California. 
We  have  also  used  seismograms  from  the  Wood-Anderson  seismograph  at  Dug- 
way,  Utah.  In  all,  471  seismograms  for  103  earthquakes  in  the  range 
0.8 %Mi  ^6.3  were  studied.  Figure  1  is  a  map  showing  the  source  regions,  sta¬ 
tion  locations,  and  travel  paths  analyzed. 

Following  Bakun  (1978),  we  have  synthesized  Wood-Anderson  seismograms 
by  transforming  the  original  time  series  into  the  frequency  domain  and  decon¬ 
volving  the  original  instrument  response,  thus  converting  the  spectrum  to 
ground  displacement.  This  complex  ground  displacement  spectrum  was  then 
convolved  with  a  Wood-Anderson  response,  and  the  resulting  spectrum 
transformed  back  to  the  time  domain  to  give  a  ’'synthethic"  Wood-Anderson 
seismogram,  i.e.  that  which  would  have  been  recorded  at  the  site  had  the  origi¬ 
nal  instrument  been  replaced  with  a  Wood-Anderson  torsion  seismograph  with 
sufficient  dynamic  range.  Maximum  zero-to-peak  amplitudes  were  determined 
from  the  Wood-Anderson  seismograms,  and  these  values  were  converted  to  a 
Ml  ( syn )  using  the  logA0  table  given  in  Richter  (1958). 

In  keeping  with  the  original  definition  of  local  magnitude,  we  have  used  epi- 
central  distance  when  determining  the  logA0  term.  Jennings  and  Kanamori 
(1983)  examined  the  effect  of  different  measures  of  distance  on  estimates  of 
near  source  Mi  and  found  that  epicentral  distance  was  not  inappropiate  when  it 
was  large  compared  to  the  fault  length.  For  all  but  the  largest  of  the  events  tnis 
assumption  is  valid  and  since  the  fault  lengths  for  the  large  earthquakes  are  not 
known,  we  used  epicentral  distance  in  those  cases  as  well. 

We  used  Mi  as  reported  by  the  Univerisity  of  California  at  Berkeley  and  the 
California  Institute  of  Technology  as  a  "true"  reference  magnitude.  When  both 
institutions  reported  a  Mi  for  a  given  earthquake  we  used  the  average  as  our 
reference,  otherwise  we  used  the  single  value  given  (which  usually  is  an  average 
over  several  stations).  Calling  the  reference  magnitude  ML(WA),  we  plotted 
Mi  (syn)  -ML  (WA  )  versus  distance.  Deviations  from  the  zero  line  in  these  plots 
indicate  differences  between  the  Great  Basin  log<40(A)  curve  and  Richter’s  (1958) 
curve  with  which  the  reference  magnitudes  were  determined.  Positive  "resi¬ 
dual”  magnitudes  reflect  larger  logA0  values  than  Richter's  while  negative  resi¬ 
duals  are  due  to  smaller  log,40  values. 
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Determination  of  logA„(A)  for  the  Great  Basin 

Figure  2a  is  the  ML  {syn )  -  Mi  { WA  )  versus  distance  plot  for  all  data  with 
Ml  S6.0.  We  observe  a  near  source  saturation  of  Mi  for  these  large  events  simi¬ 
lar  to  that  observed  for  comprable  sized  earthquakes  in  southern  California 
(Luco,  1982;  Jennings  and  Kanamori,  1983)  and  northern  Baja  California 
(Munguia  and  Brune,  1984).  Specifically,  we  find  that  for  AS  20  km  Mi  is  con¬ 
sistently  underestimated  by  a  few  tenths  to  almost  one  magnitude  unit  while  at 
greater  distances  (A  >  40  km)  we  obtain  mostly  positive  residuals. 

Munguia  and  Brune  (1984)  have  postulated  a  possible  mechanism  for  this 
observation.  They  suggest  that  low  Mi  for  large  events  recorded  at  near  dis¬ 
tances  arise  from  a  complex  near-source  spreading  out  the  waveform  in  time  as 
a  result  of  the  larger  source  size.  They  describe  an  experiment  in  which  they 
use  a  small  event  {Mi  =4.8)  as  an  emperical  Green’s  function  which  is  summed 
over  a  rectangular  fault  to  determine  the  seismogram  for  a  larger  event  at  both 
near  and  regional  distance.  They  find  that  at  greater  epicentral  distances 
(A >50  km)  the  arrivals  from  various  portions  of  the  fault  constructively  add  to 
produce  larger  amplitude  arrivals  and  hence  yield  a  greater  Wood-Anderson 
magnitude.  However,  at  near  distances  (A  <  20  km)  the  arrivals  are  not  in  phase, 
leading  to  a  lower  estimate  of  Mi .  This  same  situation  appears  to  apply  for  the 
large  events  we  have  studied. 

Figure  2b  is  a  plot  of  the  same  data  as  in  Figure  2a  except  that  Jennings 
and  Kanamori's  (1983)  modifications  to  Richter's  (1958)  attenuation  curve  have 
been  applied.  This  brings  our  data  into  agreement  with  the  California  Wood- 
Anderson  stations.  In  both  Figures  2a  and  2b,  the  data  at  the  500  to  600  km 
range  agree  reasonably  well  with  the  California  data.  This  could  lead  one  to  con¬ 
clude  that  at  those  distances  the  Great  Basin  attenuation  curve  matches  the 
southern  California  curve.  However,  we  feel  that  our  data  for  smaller  events  do 
not  support  such  a  conclusion,  as  we  discuss  below. 

In  Figure  3  we  plot  M  i  {syn)  -  M i  {WA)  versus  distance  for  all  data  with 
Mi  §5.5.  Although  there  is  considerable  scatter,  two  gross  features  are  obvious. 
First,  at  near  distances  (A<  15  km)  Richter’s  log/40  curve  results  in  an  overesti¬ 
mation  of  Mi  relative  to  the  more  distant  California  stations  by  as  much  as  one 
magnitude  unit.  Second,  at  larger  distances  (A  >  75  km)  Richter’s  curve  causes 
the  Great  Basin  estimates  of  Mi  to  be  about  0.3  magnitude  units  too  low.  We 
shall  consider  each  observation  separately. 

Figure  4a  shows  the  near  source  data  in  greater  detail  while,  for  clarity,  Fig¬ 
ure  4b  shows  the  data  for  two  individual  events.  The  reduction  in  scatter  visible 
in  Figure  4b  indicates  that  much  of  that  in  Figure  4a  is  probably  due  to  varia¬ 
tions  in  source  properties  between  different  events.  Both  plots  show  a  sys¬ 
tematic  increase  in  Mi  {syn)  -  Mi  {WA)  with  decreasing  distance,  ranging  from 
near  zero  at  12  ?  A§  30  km  up  to  about  one  at  zero  distance.  These  observations 
are  very  similar  to  those  made  by  Munguia  and  Brune  (1984)  for  comprable  sized 
earthquakes  in  southern  California  and  northern  Baja  California.  Consequently, 
we  feel  that  this  near  source  overestimation  of  Mi  is  due  to  the  fact  that 
Richter's  (1958)  attenuation  curve  is  poorly  constrained  at  close  distances  and 
not  to  any  physical  differences  between  the  western  Great  Basin  and  southern 
California. 

We  have  determined  modificatons  to  the  logA0  curve  for  A  <  20  krn  which 
remove  the  observed  bias  from  both  our  data  and  that  of  Munguia  and  Brune 
(1984,  fig.  4).  Figures  4c  and  4d  show  the  result  of  computing  Mi  with  this 
revised  curve.  We  point  out  that  the  sense  of  this  correction  for  our  close-in 
data  is  opposite  to  that  given  in  Jennings  and  Kanamori  (1983)  since  here  we  are 
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dealing  with  smaller  events  which  do  not  demonstrate  the  effect  of  finite  source 
size.  Thus,  at  near  distances  two  attenuation  curves  are  required  in  order  to 
compute  Mi  over  a  complete  magnitude  range. 

Turning  our  attention  now  to  the  data  in  the  100  to  600  km  distance  range 
(Figure  3a)  we  note  that  at  all  distances  except  at  about  600  km,  Mi  is  on  the 
average  underestimated  by  about  0.3  magnitude  units.  Those  data  at  600  km 
correspond  to  the  station  at  Dugway,  Utah  which  is  located  on  sediments.  All 
the  other  stations  at  regional  distances  are  installed  on  bedrock,  so  it  is  possible 
that  the  Dugway  "anomaly"  is  a  result  of  sediment  amplification.  We  conclude 
from  Figure  3a  that  within  the  Great  Basin  there  is  greater  attenuation  of 
seismic  energy  in  the  Wood-Anderson  passband  than  in  southern  California. 

This  result  is  not  inconsistent  with  other  geophysical  studies  in  the  Great 
Basin,  which  show  that  it  is  a  region  of  active  crustal  spreading  characterized  by 
high  heat  flow  (Sass,  et  al.,  1977)  and  a  low  velocity,  high  attenuation  upper  man¬ 
tle  (Archambeau,  et  al.,  1969).  Patton  and  Taylor  (1984)  concluded  that  the 
shear  wave  Q  was  about  100  in  the  lower  crust  and  30  in  the  upper  mantle  which 
is  less  than  the  coda  Q  of  200  Singh  and  Herrmann  (1983)  found  for  California. 
This  is  in  conflict,  however,  with  Evernden’s  (1975)  study  of  intensities  in  the 
Great  Basin  which  concluded  that  it  has  slightly  less  attenuation  than  California. 

Nevertheless,  our  data  require  a  systematic  increase  in  attenuation  relative 
to  southern  California  in  order  to  bring  our  Mi  estimates  into  agreement  with 
those  made  using  the  Berkeley  and  Pasadena  Wood-Anderson  recordings  of  the 
same  earthquakes.  We  have  done  so,  and  Figure  3b  shows  the  result  of  our 
corrections.  In  Figure  5  we  plot  our  logA0  curve  for  the  Great  Basin  together 
with  Richter's  (1958)  curve  and  Jennings  and  Kanamori’s  (1983)  modification  to 
it.  For  A  <  15  km  the  Great  Basin  curve  rises  above  Richter’s,  is  the  same  as  his 
for  15  <  A  <  35  km,  and  then  falls  below  it  by  a  constant  factor  of  0.3. 

All  of  our  work  is  predicated  on  the  assumption  that  our  reference  magni¬ 
tude  (the  average  of  the  U.  C.  Berkeley  and  Caltech  magnitudes)  is  correct,  how¬ 
ever  this  may  not  be  the  case.  Bakun  (1984a)  reports  on  a  possible  miscalibra- 
tion  of  the  Berkeley  Wood-Anderson  seismographs  in  which  it  appears  that  the 
Berkeley  instruments  operate  at  a  gain  of  2100  rather  than  2800  as  required  by 
the  definition  of  ML  .  If  this  is  the  case,  then  our  baseline  could  be  as  much  as 
0  12  magnitude  units  too  high.  This  shift  in  the  baseline  would  affect  only  the 
details  of  our  conclusions,  however.  The  near  source  overestimation  of  Mi 
observed  for  the  smaller  events  would  be  enhanced  and  although  the  underesti¬ 
mation  found  for  the  same  data  at  regional  distances  would  be  reduced,  it  would 
still  persist.  Furthermore,  since  for  most  events  we  include  the  Caltech  magni¬ 
tude  in  our  average,  the  effect  of  the  Berkeley  miscalibralion  should  be  less 
than  0.12  units. 


M0  vs  Mi  relationships  for  the  1980  Mammoth  Lakes  sequence 


Using  the  logA0(A)  curve  developed  here  for  the  smaller  events  and  Jen¬ 
nings  and  Kanamori's  (1983)  curve  for  the  larger  events,  we  recomputed  Mi  for 
71  earthquakes  of  the  1980  Mammomth  Lakes,  California  sequence.  Figure  6  is  a 
plot  of  these  revised  magnitudes  versus  the  seismic  moment  taken  from 
Archuleta,  et  al.  (1902)  and  Priestley,  et  al.  (1985).  These  moments  are  based 
on  spectral  analysis  of  the  same  digital  and  strong-moUon  recordings  from 
which  the  synthetic  Wood-Anderson  seismograms  were  derived.  Both  Archuleta, 
et  al.  (1982)  and  Priestley,  et  al.  (1985)  obtained  similar  values  of  seismic 
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moment  when  they  examined  seismograms  common  to  both  studies.  In  addi¬ 
tion,  Priestley,  et  al.  determined  the  moments  of  the  larger  (Mi^  4.0)  events 
from  regional  surface  wave  recordings  and  found  them  to  be  in  agreement  with 
those  determined  from  spectral  analysis  of  the  local  recordings.  In  the  figure, 
error  bars  denote  the  spread  in  Mi  and  Ma  for  multiple  recorded  events. 

Two  straight  lines  fit  the  data.  For  events  in  the  magnitude  range 
3.0  5  Mi  5  6.25,  the  data  are  fit  by  the  line 

log  M0  =  (1.22+0.04  )Ml  +  (16. 98±Q.  15)  (1) 

The  smaller  events  {Mi  5  1.5)  deviate  from  this  line.  For  events  in  the  magnitude 
range  0.8  5  A/),  5  3.0  the  moments  scatter  about  the  line 

logMa  =  (0.70±0.1l)Mi  +(18.263:0.24) 

however  it  is  strongly  affected  by  the  four  smallest  events.  There  are  very  few 
recordings  for  these  four  events  and  it  could  be  easily  argued  that  neither  their 
moments  nor  their  magnitudes  are  well  determined. 

The  M0  versus  Mi  relationship  given  in  equation  (l)  can  be  compared  to 
that  given  by  Archuleta,  et  al.  (1982)  which  used  data  primarily  in  the  range 
3.0  Is  Mi  5  5.0  with  only  one  event  above  Mi  =5.0.  Our  data,  on  the  other  hand, 
span  the  range  0.8  5  Mi  5  6.3  with  eight  events  above  Mi  =5.0.  Archuleta,  et  al. 
(1982)  found  that  the  M0  versus  Mi  relationship  for  the  40  events  for  which  they 
had  Berkeley  Mi  determinations  was,  for  the  range  2.9%  Mi  5  6.2, 

log  M0  =  (0.96+0.06)^  +  (18. 14±0.23) 

but  for  the  larger  events  with  more  reliable  Mi  estimates  (3.5  5  ML  5  6.2)  it  was 

log M Q  =  (1.05±0.08 )Ml  +  ( 17.  76±0.33) 

The  greater  slope  in  equation  (l)  is  due  to  the  larger  moments  which  we  used  for 
events  with  A//,  >5.0.  For  example,  for  the  May  27,  1980  event  at  14:50, 
Archuleta,  et  al.  (1982)  found  a  seismic  moment  of  2.33x’024  dyne-cm  based  on 
two  recordings  whereas  Priestley,  et  al.  (1985)  obtainec  a  value  of  7.24xlOa4 
dyne-cm  when  using  nine  recordings. 

In  comparing  our  data  with  other  regions  of  California,  the  slope  of  our 
curve  is  less  than  the  least  squares  fit  Archuleta,  et  al.  (1982)  made  to  the  data 
in  Thatcher  and  Hanks  (1973),  but  is  essentially  identical  to  the  relation  found 
by  Bakun  and  Lindh  (1977) 

log  M0  =  (1.21  +  0.03  )Ml  +  (17.02+0,07) 

for  the  Oroville,  California  earthquake  sequence.  Oroville  is  located  in  the 
western  Sierra  Nevada  in  a  tectonic  environment  which  is  more  similar  to  that  of 
Mammoth  Lakes  than  the  southern  California  region  studied  by  Thatcher  and 
Hanks  (1973). 

In  a  recent  study  of  five  diffeient  source  regions  (Parkfield,  San  Juan 
Bautista,  the  Sargent  fault,  Coyote  Lake,  and  the  Livermore  Valley),  Bakun 
(1984b)  detected  an  inflection  in  the  log  M0  versus  Mi  plot  at  about  Mi  =3.  For 
the  events  in  the  range  1.5  5  Mi  5  3.5  he  found  the  relationship  to  be 

log  M0  =  1.2  Ml  +  17 

which  is  essentially  the  same  relationship  given  in  equation  (l)  based  on  data  in 
the  range  3,05  ML  5  6.25.  However,  for  events  in  the  similar  range  3.5  5  Mi  5  6.25 
Bakun  (1984b)  found  the  relationship  to  be 

log  M0  =  l.bMi  +  16 
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Recently,  Hanks  and  Boore  (1984)  summarized  a  number  of  M0  versus  Mi 
relationships  published  for  central  and  southern  California  and  compared  them 
with  model  calculations  based  on  Boore  (1983).  In  Boore  (1983),  seismograms  of 
hypothetical  earthquakes  are  generated  by  forcing  the  amplitude  spectrum  of 
band-limited,  white  Gaussian  noise  to  be  that  predicted  by  the  Brune  (1970, 
1971)  source  model.  Boore  assumed  a  constant  stress  drop  (of  100  bars),  thus 
restricting  the  spectrum  scaling  law  to  be  a  function  of  seismic  moment  alone. 

Hanks  and  Boore  (1984)  computed  M0  and  Mi  from  these  "seismograms" 
and  found  that  predicted  values  of  M0  versus  Mi  fit  their  data  up  to  about 
Mi  =  6.0  at  which  point  their  Mi  values  exceeded  those  of  the  model.  In  Figure  6 
we  have  plotted  the  same  theoretical  values  as  given  in  Hanks  and  Boore  and 
find  that  the  predictions  agree  with  our  data  for  small  magnitude  events  (if  we 
neglect  the  four  smallest  earthquakes  mentioned  earlier)  and  depart  from  our 
observations  at  around  Mi  =  5.0.  Hanks  and  Boore  (1984)  used  teleseismic  esti¬ 
mates  of  seismic  moment  for  the  four  largest  Mammoth  Lakes  events  which  are 
20  to  30  percent  greater  than  the  estimates  Priestley  et  at.  (1985)  obtained 
through  spectral  analysis  of  local  recordings  of  these  same  events.  Had  we  used 
the  larger  moments,  the  predicted  M0  versus  Mi  values  would  deviate  less  from 
the  observations.  The  fact  that  Boore's  (1983)  model  (which  assumes  constant 
stress  drop)  and  our  data  do  not  agree  for  large  Mi  implies  that  stress  drop  for 
the  Mammoth  Lakes  events  is  not  constant  but  rather  increases  with  magnitude. 
This  is  in  agreement  with  results  from  spectral  analysis  (Archuleta,  et  at.  1982; 
Priestley,  et  al..  1985)  which  indicate  that  the  largest  events  had  stress  drops  on 
the  order  of  1000  bars. 


Summary 

Comparison  of  peak  Wood-Anderson  response  measurements  in  the  Great 
Basin  to  those  in  California  for  the  same  events  indicate  that  there  is  a  near 
source  Mi  saturation  for  events  with  Mi  §  6.0.  For  smaller  events  {Mi  g  5.5)  local 
magnitude  is  increasingly  overestimated  at  distances  less  than  about  15  km, 
probably  due  to  the  poor  constraint  on  Richter’s  (1958)  attenuation  curve.  At 
larger  distances,  Mi  in  the  Great  Basin  is  underestimated  by  approximately  0.3 
units,  suggesting  greater  attenuation  than  in  California.  We  have  modified 
Richter’s  (1958)  attenuation  curve  so  as  to  remove  the  above  bias  for  Great 
Basin  seismograms.  A  revised  moment-magnitude  relationship  for  the  1980 
Mammoth  Lakes  earthquakes  based  on  the  modified  curve  is  similar  to  that 
found  for  Sierra  Nevada  earthquakes.  The  relationship  obtained  suggests  a  con¬ 
stant  stress  drop  of  100  bars  for  smaller  Mammoth  earthquakes  but  increasingly 
larger  stress  drops  for  events  above  Mi  =  5. 
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Figure  captious 


Figure  1.  Map  of  the  western  United  States  showing  the  Mammoth  Lakes,  Califor¬ 
nia  source  region  (enlargement,  epicenters  plotted  as  dots),  seismic  sta¬ 
tions,  and  ray  paths  for  the  data  used  in  this  study.  The  heavy  line  del¬ 
ineates  the  border  of  the  Basin  and  Range.  The  triangles  indicate  seismic 
stations  used  to  determine  the  Great  Basin  attenuation  curve,  the  two  solid 
squares  give  the  locations  of  Berkeley  and  Pasadena,  California. 

Figure  2.  Mi  ( syn  )  -  Mi  (WA)  versus  distance  for  all  events  with  ML  (WA  )  2  6.0.  a) 
Mi  (syn)  computed  using  Richter's  (1958)  log/40  values,  b)  Mi(syn)  com¬ 
puted  using  Jennings  and  Kanamori’s  (1983)  modifications  to  Richter’s 
values. 

Figure  3  Mi  (syn)  -  Mi  (WA  )  versus  distance  for  all  events  with  Mi  (iF4)^5.5,  a) 
Mi  (syn)  computed  using  Richter’s  (1958)  log.40  values,  b)  Mi  (syn)  com¬ 
puted  using  the  Great  Basin  log.40  curve  developed  in  this  study. 

Figure  4.  a)  The  same  data  as  in  Figure  3a  with  epicentral  distance  under  30 
km.  b)  Two  individual  events  from  the  data  set  shown  in  Figure  4a.  Event 
dates  and  origin  time  are  given  In  the  upper  right  corner,  c  and  d)  The 
same  data  as  in  (a)  and  (6)  except  with  Mi  (syn)  computed  using  the  Great 
Basin  log Aa  curve. 

Figure  5.  Comparison  of  the  log^0  curve  obtained  for  the  Great  Basin  with 
Richter’s  (1958)  southern  California  curve  and  Jennings  and  Kanamori’s 
(1983)  modification  to  it. 

Figure  8.  Log  Ma  versus  Mi  for  earthquakes  in  the  1980  Mammoth  Lakes,  Califor¬ 
nia  sequence.  Stars  indicate  theoretical  values  given  in  Hanks  and  Boore 
(1984)  assuming  a  constant  stress  drop  of  100  bars. 
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